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Abstract 

This  research  is  a  foundational  study  of  conformal,  load-bearing  antenna  arrays 
embedded  into  the  wing  structure  of  a  joined-wing  aircraft.  It  is  a  multidisciplinary 
effort  that  touches  on  the  aerodynamic,  structural,  and  electromagnetic  design  con¬ 
siderations  that  stem  from  this  unique  type  of  sensor  integration.  The  antenna 
performance,  Finite  Element  Model  (FEM),  and  control  surface  effectiveness  are 
investigated.  The  theory  describing  an  ensemble  of  dipole  antenna  elements  that 
conform  to  the  shape  of  a  section  of  the  joined  wing  is  developed.  The  far  field,  free 
space  radiation  pattern  of  the  sensor  is  then  analyzed  for  a  wing  that  is  deflected  due 
to  typical  aerodynamic  loading.  This  pattern  is  compared  to  the  same  antenna  when 
the  wing  is  not  deformed.  A  FEM  of  the  antenna  elements  is  created  and  incorpo¬ 
rated  into  the  full  FEM  of  the  joined-wing  aircraft  allowing  its  structural  impact  to 
be  realized.  Based  on  the  positioning  of  these  large  sensor  arrays,  control  surfaces 
are  placed  and  examined  to  achieve  the  proper  handling  capabilities  necessary  for 
this  type  of  aircraft. 

The  results  of  this  study  show  that  wing  deflections  due  to  typical  aerodynamic 
loads  produce  significant  disturbances  to  the  radiation  pattern  of  a  conformal  an¬ 
tenna  when  end-fire  phasing  is  applied.  Active  compensation  for  wing  deformation 
will  have  to  be  applied  to  correct  the  beam  steering.  The  basic  array  theory  provides 
a  useful  tool  for  analyzing  the  radiation  pattern  of  an  array  that  conforms  to  the  sur¬ 
face  of  an  undeformed  wing.  On  wings  deformed  due  to  typical  aerodynamic  loads, 
the  main  beam  elevation  pointing  angle  varies  considerably  between  the  basic  array 
theory  and  NEC-Win  Plus+™  solutions.  Also,  the  reversal  speed  of  an  outboard 
aileron  is  determined. 
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STRUCTURALLY  INTEGRATED  ANTENNAS  ON  A 


JOINED-WING  AIRCRAFT 


I.  Introduction 

This  research  effort  focused  on  load-bearing,  conformal  antennas  embedded 
into  the  skin  of  a  joined-wing  type  of  aircraft.  It  was  part  of  the  ongoing  work 
described  in  [2],  One  of  the  main  areas  of  interest  was  determining  the  pointing 
accuracy  of  the  antenna  array  as  the  wing  was  deflected  due  to  typical  aerodynamic 
loads  acting  on  the  lifting  surfaces. 

A  simplified,  composite  FEM  of  Northrop  Grumman’s  Sensorcraft  -  Conformal 
Load  Bearing  Antenna  Structure  [15]  (S-CLAS)  was  created.  This  new  element  was 
then  substituted  into  the  joined-wing  model  in  place  of  the  elements  on  each  wing 
surface  where  the  antenna  arrays  will  be  located. 

Placement  of  the  antennas  is  an  important  issue  when  designing  the  aircraft 
due  to  the  fact  that  the  conformal  antennas  utilize  the  entire  chord  of  the  wing. 
Therefore,  control  surfaces  cannot  be  placed  in  any  spanwise  section  of  the  wings 
where  the  antennas  arc  located.  The  size  of  the  antennas  is  performance  driven,  but 
there  are  size  limitations  based  on  the  amount  of  control  surfaces  needed  to  effectively 
turn  and  trim  the  joined-wing  aircraft.  Thus,  the  control  surface  effectiveness  is 
another  area  of  interest. 

1.1  Background 

Unmanned  Arial  Vehicles  (UAVs)  have  become  an  integral  part  of  the  aero¬ 
nautical  engineer’s  vocabulary.  These  aircraft  give  the  ability  to  maintain  constant 
surveillance  of  ground  activities.  The  notional  sensorcraft  is  a  UAV  concept  in  which 

the  Air  Force  Research  Laboratory  (AFRL)  is  fusing  many  of  its  state-of-the-art 
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technologies  into  a  single  aircraft  [2],  The  AFRL  concept  uses  a  joined  front  and  aft 
wing  and  conformal  antennas  integrated  into  the  wing  structure.  The  wings  connect 
in  a  diamond  pattern  when  viewed  from  the  front  and  from  above  as  shown  in  Fig¬ 
ure  1.1  [2],  This  diamond  shape  allows  the  sensorcraft  to  view  360°  of  its  environment 
in  the  azimuth  plane,  thus  allowing  fewer  potential  targets  to  go  undetected. 


REAR  (right  wing  up) 


SIDE 

Figure  1.1  Joined-Wing  Configuration 

Table  1.1  [2]  shows  the  geometric  parameters  of  the  current  configuration.  One 
can  see  that  this  aircraft  has  a  very  long  tip-to-tip  wingspan  of  approximately  64.5 
meters  (212  ft).  A  span  of  this  magnitude  makes  the  structure  quite  flexible,  and 
this  is  a  concern  due  to  the  embedded  antennas.  This  type  of  design  is  not  optimal 
when  you  are  trying  to  keep  a  radar  beam  focused  on  a  target.  The  amount  of  wing 
flexing  allowed  by  the  sensors  tends  to  drive  the  structure  to  be  stiffen  The  sensors 
are  discussed  in  more  detail  later  in  this  section.  An  overview  of  the  structural 
aspects  of  the  joined-wing  is  first  presented. 

Many  joined-wing  designs  have  been  studied  over  the  past  17  years.  The  first 
joined-wing  study  was  presented  at  a  conference  by  J.  Wolkovitch  [25]  in  1985.  Livne 
presented  an  extensive  survey  of  papers  stemming  from  Wolkovitch’s  work  and  oth- 

1-2 


Table  1.1  Baseline  Configuration  Parameters 


DESCRIPTION 

LABEL 

VALUE 

Inboard  Span 

sib 

26.00  m 

Outboard  Span 

Sob 

6.25  m 

Forward  Root  Chord 

crf 

2.50  m 

Aft  Root  Chord 

Cra 

2.50  m 

Mid-Chord 

Cm 

2.50  m 

Tip  Chord 

ct 

2.50  m 

Forward-  aft-x-off  set 

Xfa 

22.00  m 

Forward- aft-z-offset 

Zfa 

7.00  m 

Inboard  Sweep 

A  ib 

30  deg 

Outboard  Sweep 

A ob 

30  deg 

Airfoil 

LRN-1015 

Calculated  Planform  Area 

S 

145.0  m2 

Calculated  Wing  Volume 

52.2  m3 

ers  up  to  the  year  2001  [7].  This  survey  is  a  valuable  document  for  anyone  doing 
research  in  any  area  of  joined- wing  design.  One  pertinent  statement  in  his  paper  was 
the  fact  that  antenna/lifting  surface  interactions  should  be  studied.  No  past  work 
accomplished  in  this  area  was  referenced.  G.  Reich  et.  al.  presented  their  work  [17], 
which  minimized  deformation  in  the  antenna  by  utilizing  Active  Aeroelastic  Wing 
technology,  but  evaluating  the  antenna  performance  under  these  conditions  was  not 
the  focus  of  their  research.  They  did  not  model  the  antennas  in  this  study  as  load- 
bearing  parts  of  the  structure,  and  therefore  no  stiffness  properties  were  included  for 
the  antennas. 

A  conformal  array  type  of  antenna  is  necessary  to  allow  the  antennas  to  be 
embedded  into  the  wing  structure  while  maintaining  a  thin,  efficient  wing.  The  array 
then  becomes  a  load-bearing  member.  In  the  present  study  the  array  is  used  in  an 
end-fire  configuration,  with  the  main  beam  pointing  normal  to  the  leading  edge  of 
the  wing.  Proper  phasing  of  the  antenna  elements  in  the  array  allows  the  beam  to 
be  steered  in  a  plane  extending  off  the  end  of  the  array.  R.  Mailloux  [9]  explains 
that  conformal  arrays  may  be  treated  as  planar  arrays  if  the  length  of  the  array  is 
much  less  than  the  radius  of  curvature  of  the  surface.  Although  this  would  make 
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the  problem  much  simpler  and  might  be  relevant  for  a  static  scenario,  the  arrays 
in  this  case  should  be  modeled  as  non-planar  because  the  interest  is  in  determining 
the  disturbance  caused  by  wing  deformations.  In  this  dynamic  environment  the 
orientation  of  the  array  elements  will  be  constantly  changing  due  to  aerodynamic 
effects  on  the  flexible  wing,  thus  changing  the  pointing  direction  of  each  element 
in  the  array.  Therefore,  modeling  the  antenna  as  a  conformal  array  right  from  the 
beginning  is  the  best  approach.  Some  simplifying  assumptions  about  the  interactions 
between  elements  and  the  array  geometry  are  necessary  in  order  to  maintain  the 
scope  of  this  research  and  adapt  planar  array  theory  concepts  [1,21]  to  this  particular 
conformal  array.  These  tradeoffs  must  be  examined  so  that  a  solution  is  obtained 
that  is  a  realistic  approximation. 

Much  of  the  textbook  literature  about  arrays  [1,  9,  21]  focuses  primarily  on 
linear  arrays.  The  basic  concepts  from  linear  array  theory  are  then  used  to  construct 
planar  arrays.  W.  Kummer  [5]  comments  that  a  planar  array  is  essentially  a  string 
of  linear  arrays  placed  next  to  each  other  forming  a  two  dimensional  array.  Due  to 
the  complexity  of  conformal  arrays,  many  of  the  current  textbooks  focus  on  confor¬ 
mal  arrays  that  span  the  circumference  of  a  cylinder.  Although  the  geometry  of  a 
cylindrical  array  is  somewhat  simple,  even  this  case  requires  complex  phasing  of  each 
element  in  the  array  in  order  to  get  proper  beam  shaping.  When  one  transitions  to 
conformal  arrays,  where  each  element  may  be  pointing  in  a  different  direction,  beam 
forming  and  beam  steering  become  a  complicated  phasing  issue. 

When  defining  the  antenna,  there  are  typically  two  approaches.  If  the  dimen¬ 
sions  of  the  array  (i.e.  the  number  of  elements  and  elemental  spacing  or  the  length  of 
the  array)  are  set  as  requirements,  then  the  beamwidth  and  directivity  are  derived. 
Otherwise,  the  requirements  are  placed  on  the  beamwidth  and  directivity.  Balanis’ 
textbook  [1]  gives  a  graphical  representation  of  how  the  half-power  beamwidth  (a 
typical  measurement  for  describing  beamwidth)  varies  with  the  length  of  the  array, 
which  is  measured  in  terms  of  wavelengths.  In  the  initial  effort  of  the  research  the 
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dimensions  of  the  array  were  specified  at  approximately  five  feet  in  the  chordwise 
direction  and  42  feet  spanwise. 

1.2  Problem  Statement 

A  conformal  load-bearing  antenna  is  a  new  concept  in  the  sensor  community. 
This  type  of  antenna  replaces  the  skin  material  on  portions  of  the  wings  and  therefore 
carries  the  aerodynamic  loads  that  the  skin  would  normally  experience.  As  the  wings 
deform  due  to  typical  flight  conditions,  the  antenna  is  also  bending  and  twisting 
away  from  its  initial  configuration.  There  are  currently  no  methods  that  quantify 
the  effects  of  this  deformation  on  the  conformal  antenna  performance,  such  as  the 
changes  in  the  main  acquisition  beam’s  magnitude  and  direction.  A  method  should 
be  developed  that  quickly  and  accurately  measures  the  changes  in  these  parameters. 

1.3  Research  Objectives 

The  goal  of  this  research  was  to  quantify  the  effect  that  wing  deformation 
has  on  the  pointing  accuracy  of  conformal  load-bearing  antennas  on  a  joined-wing 
aircraft.  Three  objectives  were  accomplished: 

1.  A  finite  element  model  of  the  conformal  was  developed  and  incorporated  into 
the  full  aircraft  FEM. 

2.  A  theoretical  antenna  modeling  approach  for  this  antenna  configuration  was 
developed  to  quantify  the  change  in  antenna  performance  due  to  deformations. 
This  approach  was  verified  using  a  commercially  available  Method  of  Moments 
software  package. 

3.  The  control  surface  effectiveness  was  determined  given  the  size  and  location  of 
the  antenna  arrays. 
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1.4  Methodology 

In  an  effort  to  control  the  scope  of  this  research,  an  array  of  horizontally 
polarized,  half-wavelength  dipoles  were  used  in  the  antenna  analysis  to  model  the 
actual  conformal  load-bearing  antenna  element.  This  allowed  the  basic  theory  of  the 
antenna  ensemble  to  be  developed  and  verified.  A  MATLAB  [10]  code  was  developed 
to  plot  the  radiation  pattern  of  the  array  based  on  the  location  and  orientation  of 
each  antenna  element.  A  commercial  software  package,  NEC-Win  Plus+™  [14],  was 
then  used  to  verify  the  analytical  model.  This  software  package  utilizes  the  Method 
of  Moments  solution  technique  to  analyze  the  antenna  and  determine  its  radiation 
pattern.  Figure  1.2,  created  in  NEC-Win  Plus+™,  shows  one  chordwise  array  of 
dipoles  conforming  to  the  shape  of  the  airfoil. 


Leading  Edge  Dipole 


Trrailing  Edge  Dipole 


Figure  1.2  One  chordwise  array  of  dipoles 

The  proposed  antennas  for  the  sensorcraft  are  conformal  arrays  that  are  inte¬ 
grated  into  the  skin  of  the  fore  and  aft  wings.  The  antennas  and  inboard  control 
surfaces  are  conhned  in  the  spanwise  direction  to  fit  between  the  fuselage  and  the 
wing  joint.  The  amount  of  wing  space  available  for  the  antennas  will  drive  the  point¬ 
ing  accuracy  of  the  main  beam,  the  maximum  beam  coverage  area,  and  ultimately, 
the  ability  to  find  and  track  targets. 

A  simple  model  of  the  sensors  was  generated  to  incorporate  the  different  layers 
of  materials  that  make  up  the  antenna’s  structure.  This  allowed  a  more  accurate 
model  of  the  wing  structure  to  be  generated.  This  also  affects  the  overall  stiffness 
of  the  structure,  which  in  turn  affects  the  amount  of  deflection  in  the  wing.  The 
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simplified  model  consisted  of  an  electromagnetically  transparent  material  called  As- 
troquartz,  honeycomb  core  material,  and  graphite  epoxy,  which  acts  as  the  main 
load-bearing  component  of  the  antenna  element.  This  element  model  replaces  the 
current  aluminum  elements  in  the  skin  of  the  wings  where  the  antenna  array  is  to 
be  situated.  The  representation  of  the  FEM  is  shown  in  Figure  1.3. 


Astroquartz _ 


Honeycomb 
Core  Structure 


Graphite  Epoxy 


Figure  1.3  Simplified  finite  element  model  of  antenna 


A  baseline  configuration  of  the  array  was  generated  to  correspond  to  the  an¬ 
tenna  orientation  assuming  the  wing  was  not  deformed.  The  antenna  modelling  soft¬ 
ware  then  generated  the  beam  pattern.  Phasing  of  the  array  elements  was  performed 
to  generate  the  desired  baseline  beam  orientation.  Once  this  baseline  was  established, 
specific  load  configurations  were  applied  to  the  lifting  surfaces.  Wing  deformations 
were  then  be  generated  on  a  fully  stressed  design  [2]  by  an  integrated  software  en¬ 
vironment  utilizing  the  Adaptive  Modeling  Language  (AML),  MSC.NASTRAN  and 
PanAir  [10,  16,  22],  These  deformations  were  used  to  determine  the  new  location 
and  slope  of  each  element  of  the  array.  Figure  1.4  shows  the  typical  deformed  shape 
of  the  array  compared  to  the  undeformed  shape.  A  new  beam  pattern  was  then 
generated.  This  process  was  repeated  for  the  various  load  configurations.  The  re¬ 
sults  were  compared  to  the  baseline  results  to  determine  the  beam  pointing  errors 
produced  by  the  wing  deformations. 

Based  on  the  location  of  the  antenna  arrays  along  the  span  of  the  wing  sections, 
the  placement  and  effectiveness  of  the  control  surfaces  varies.  The  effectiveness  was 
checked  with  just  one  control  surface  outboard  of  the  wing  joint.  The  process  was  not 
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Figure  1.4  Deformed  (upper)  and  undeformed  (lower)  configurations  of  the  an¬ 
tenna  array 

optimized  to  determine  the  best  configuration  of  the  antennas  and  control  surfaces. 
The  antenna  requirements  drive  the  size  of  the  array.  The  arrays  must  be  placed 
inboard  of  the  joint  on  both  the  front  and  aft  wings.  It  was  desired  to  have  maximum 
control  surface  effectiveness  while  maintaining  a  360°  beam  coverage  area  around  the 
perimeter  of  the  aircraft. 

1.5  Assumptions/Limitations 

There  were  many  aspects  of  this  research  that  could  have  been  explored  in 
great  depth,  but  this  study  was  meant  to  set  the  foundation  and  basic  theory  of 
the  interactions  between  the  antennas  and  aircraft  structure.  Therefore,  several 
assumptions  were  made  in  order  to  control  the  scope  of  this  effort.  These  assumptions 
are  broken  down  into  three  basic  categories  based  on  the  outline  of  research  objectives 
stated  earlier:  structural,  antenna,  and  control  surface  effectiveness. 

1.5.1  Structural.  A  simplified  finite  element  model  of  the  S-CLAS  incor¬ 
porating  only  the  main  load  bearing  layers  of  each  element  in  the  array  was  used. 
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By  making  this  simplification,  the  number  of  layers  of  material  in  each  element  de¬ 
creased  from  approximately  20  in  the  actual  antenna  down  to  only  five  in  the  FEM. 
This  lay-up  is  shown  in  Figure  1.3  in  the  previous  section. 

The  top  layer  of  the  FEM  was  assumed  to  be  made  of  Astroquartz  II  RS12-B, 
which  acts  as  an  electromagnetically  transparent  layer  to  allow  most  of  the  antenna’s 
energy  to  pass  through  unaffected.  The  second  and  fourth  layers  down  from  the 
top  are  made  of  honeycomb  core  material.  These  layers  act  essentially  as  spacers 
between  the  various  elements  of  the  antenna,  but  they  also  carry  some  shear  load. 
The  third  and  fifth  layers  are  the  main  load-bearing  parts  of  the  antenna.  They 
are  two  composite  graphite  epoxy  (IM7/977-3)  layers  that  act  as  a  backplane  (third 
layer)  and  a  groundplane  (bottom  layer)  in  the  antenna  structure.  Virtually  the 
entire  load  is  carried  in  these  two  layers.  These  properties  for  these  materials  are 
tabulated  in  Chapter  3. 

1.5.2  Antenna.  There  were  several  assumptions  made  with  respect  to  the 
antenna  model.  First,  the  actual  conformal  antenna  element  was  much  too  compli¬ 
cated  of  a  structure  to  accurately  model  in  detail  and  incorporate  into  this  study, 
given  the  time  constraints  and  research  objectives.  Therefore  horizontal  dipoles  were 
used  instead  of  the  actual  elements.  This  allowed  the  theory  to  be  more  easily  de¬ 
veloped  and  verified.  These  dipoles  gave  the  same  polarization  expected  from  the 
S-CLAS  [15].  The  baseline  configuration  of  the  antenna  was  assumed  to  be  a  linear 
translation  out  the  span  of  one  chordwise  array  of  dipoles.  The  dipoles  along  one 
chord  were  positioned  according  to  the  geometric  center  of  the  actual  square  ele¬ 
ments  based  on  the  assumed  airfoil,  LRN-1015  (see  Appendix  A).  Their  orientation 
was  achieved  assuming  a  simple  translation  in  the  chordwise  direction  due  to  wing 
twist  and  a  central  difference  method  in  the  spanwise  direction  to  account  for  the 
rotation  of  the  elements  due  to  the  wing  deflection.  Each  chordwise  array  was  as¬ 
sumed  to  be  rotated  as  a  single  unit  equal  to  the  amount  of  rotation  in  the  center 
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element  of  that  particular  chordwise  array.  It  was  also  assumed  that  each  individual 
dipole  remained  linear  and  was  not  allowed  to  bend. 

When  developing  the  equation  of  the  magnitude  of  the  electromagnetic  field 
(E-field),  far  held  assumptions  were  used.  This  means  that  the  vectors  pointing 
from  the  center  of  each  element  to  the  far  held  observation  point  were  assumed  to  be 
parallel.  Also,  it  was  assumed  that  the  magnitude  and  phase  of  the  E-held  radiated 
from  each  element  could  be  summed  in  the  far  held.  Therefore  no  mutual  coupling 
between  elements  was  accounted  for  in  this  approximation. 

Another  limitation  was  that  this  study  did  not  take  into  account  that  the 
surrounding  wing  structure  would  have  an  effect  on  the  antenna  performance.  For 
the  purposes  of  this  study,  it  was  sufficient  to  place  the  elements  in  the  correct 
geometrical  configuration  according  to  the  position  and  geometry  of  the  wing,  but 
assume  that  they  were  operating  in  free-space  with  no  degradation  caused  by  the 
surrounding  structure. 

The  current  distribution  along  each  dipole  was  assumed  to  be  a  cosine  taper 
with  zero  current  at  both  ends  and  a  maximum  value  of  one  at  the  center.  This  cur¬ 
rent  is  more  accurate  than  a  constant  current  because  the  current  on  a  real  element 
must  go  to  zero  at  the  ends.  The  continuous  curvature  is  also  more  practical  than 
a  triangularly  shaped  current.  A  closed  form  solution  was  then  used  for  the  Fourier 
Transform  of  this  current.  An  explanation  of  the  Fourier  Transform  is  outlined  in 
Appendix  B. 

1.5.3  Control  Surface  Effectiveness.  The  control  surface  effectiveness  was 
derived  by  using  a  program  called  M  SC.  Flight  Loads  [11]  and  NASTRAN.  Flight- 
Loads  was  used  to  create  an  aerodynamic  model  with  control  surfaces  and  spline 
this  together  with  the  structural  model.  The  control  surface  was  assumed  to  be 
the  rear  half  of  the  outboard  portion  of  the  joined  wing.  The  aero/structure  model 
was  then  processed  using  NASTRAN  to  derive  the  values  necessary  for  calculating 
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the  control  surface  effectiveness.  FlightLoads  utilizes  a  method  compatible  with  the 
Doublet-Lattice  method  for  subsonic  flows  and  ZONA51  for  supersonic  flows.  These 
methods  represent  the  lifting  surfaces  as  flat  panels  that  are  typically  aligned  with 
the  flow.  For  further  information  on  these  calculations,  the  reader  is  referred  to 
Reference  [11]. 

1.6  Implications 

This  research  allows  for  a  quick  approximation  of  the  E-field  radiated  from  an 
antenna  ensemble  of  dipole  elements  that  conform  geometrically  to  a  section  of  the 
wing  surface  given  the  coordinates  of  the  center  points  of  the  elements.  Plots  are 
generated  for  the  E-field  magnitude  for  a  baseline,  undeformed  wing  and  a  deformed 
wing,  and  the  percent  change  in  the  maximum  amplitude  is  calculated.  This  gives  the 
user  an  indication  of  the  effect  that  wing  deformation  has  on  the  pointing  accuracy 
of  the  given  antenna  ensemble. 
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II.  Literature  Review 


2.1  The  Beginnings  of  the  Joined-  Wing  Aircraft 

Wolkovitch  first  developed  the  joined-wing  aircraft  in  1976  [23].  The  basic  idea 
of  this  design  was  that  the  aft  wing  was  swept  forward  to  connect  with  the  trailing 
edge  of  the  front  wing.  This  configuration  creates  a  diamond  shape  when  viewing 
the  aircraft  from  a  front  or  top-down  perspective.  Wolkovitch  also  presented  similar 
designs  in  1982  [24],  and  in  1986  [25]  he  compared  the  performance  of  a  joined  wing 
to  a  conventional  wing  on  a  transport  aircraft  .  He  claimed  that  a  joined  wing  had 
several  advantages  over  conventional  wings,  namely  low  induced  drag,  high  trimmed 
maximum  coefficient  of  lift,  and  good  stability  and  control. 

Wolkovitch  also  found  three  other  artifacts  that  are  relevant  to  this  study. 
First,  he  noted  that  the  out-of-plane  component  of  lift  acting  on  the  wing  structure 
bends  the  wing  about  an  axis  that  is  tilted  forward  (see  Figure  2.1).  This  phe¬ 
nomenon  puts  a  downward  pitching  moment  on  the  front  wing,  inducing  a  twist 
that  is  not  present  with  conventional  aircraft.  This  twisting  motion  is  an  important 
consideration  for  antenna  arrays  that  conform  to  the  surface  of  the  wing  because  it 
adds  another  element  of  wing  deformation  besides  the  typical  bending  motion. 

Secondly,  Wolkovitch  found  from  wind  tunnel  tests  that  several  configurations 
of  control  surfaces  were  effective  for  pitch  attitude  control.  He  found,  though,  that 
using  rear  wing  elevators  does  not  create  much  additional  lift  on  the  front  wing.  This 
is  an  important  consideration  for  the  sensorcraft  concept  because  wing  “real  estate” 
is  limited  due  to  the  large  size  and  number  of  sensors.  He  also  noted  that  ailerons 
were  more  effective  on  the  rear  wing,  especially  if  the  front  wing  had  a  considerable 
dihedral  angle. 
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Figure  2.1  Modified  wing-box  structure  to  compensate  for  tilted  bending  axis  of  a 
joined  wing  showing  components  of  lift  and  induced  pitching  moment. 

2.2  Developing  the  Concept 

Smith,  Cliff,  and  Kroo,  in  conjunction  with  NASA  Ames  Research  Center, 
designed  a  flight  demonstrator  aircraft  with  joined  wings.  Their  objectives  were 
to  demonstrate  good  handling  characteristics  and  validate  the  existing  joined  wing 
design  methods.  Their  JW-1  design,  a  modified  NASA  AD-1  aircraft,  was  based  on 
the  idea  that  the  joined-wing  aircraft  would  be  most  beneficial  as  a  medium-range 
transport  due  the  potential  weight  savings  offered  by  this  type  of  configuration.  Due 
to  the  requirement  for  good  handling  qualities,  an  emphasis  was  placed  on  improving 
the  stall  characteristics  of  the  aircraft.  To  achieve  this  goal,  Smith,  Cliff,  and  Kroo 
decided  to  compromise  cruise  performance  and  add  twist  to  the  inboard  section  of 
the  front  wing.  A  small  amount  of  additional  twist  was  then  needed  in  the  aft  wing 
in  order  to  maintain  an  adequate  amount  of  lift.  By  adding  the  twist  in  the  wing 
to  ensure  correct  stall  characteristics,  an  11%  penalty  was  paid  in  the  amount  of 
induced  drag  during  cruise  [19].  Vortilons  were  also  added  to  the  configuration  to 
counter  an  upward  pitch  that  occurred  at  stall  conditions.  The  vortilons  remained 
on  the  aircraft  for  the  remainder  of  the  wind  tunnel  tests,  including  those  for  control 
surface  effectiveness. 
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Various  control  surface  configurations  were  investigated  in  this  study.  Control 
surfaces  were  segmented  along  the  trailing  edge  of  both  the  front  and  aft  wings 
to  allow  each  segment  to  act  either  independently  or  in  concert.  For  example,  by 
linking  the  two  inboard  segments  (one  on  the  front  wing  and  one  on  the  aft  wing) 
and  moving  them  equally  in  opposite  directions,  they  found  that  pitch  control  could 
be  achieved. 

A  one-sixth  scale  wind  tunnel  model  of  the  JW-1  design  was  built 
and  tested  in  the  NASA  Ames  12-ft  wind  tunnel.  Wind-tunnel-test  data 
were  used  to  assess  performance,  stability,  and  control,  and  to  confirm 
that  the  relatively  simple  design  tools  used  for  the  JW-1  are  suitable  for 
designing  joined  wings.  [19] 

A  study  was  also  performed  [6]  using  the  JW-1  in  1990,  by  Lin,  Jhou,  and 
Stearman,  to  determine  the  amount  of  influence  that  the  type  of  joint  has  on  the 
strength  and  stiffness  of  a  joined  wing.  Eight  different  joint  configurations  were  in¬ 
vestigated  in  this  study  including  a  ball  joint,  pin  joints,  universal  joint,  link  joints, 
and  a  rigid  joint.  They  were  interested  in  the  influence  of  the  joint  on  the  distri¬ 
bution  of  stress  and  aeroelastic  stability.  All  but  the  link  joints  and  the  telescoping 
joint  reduced  the  bending  moment  in  the  front  wing.  “Only  the  rigid  joint  satisfies 
all  of  the  other  criteria  of  the  minimum  strain  energy,  most  uniform  strain  energy 
density,  and  minimum  level  of  joint  reaction  loads  to  accomplish  the  roll  type  bend¬ 
ing  moment  reduction.”  [6]  This  study  was  somewhat  limited  due  to  the  fact  that 
linear  modeling  was  employed  and  only  under  cruise  flight  conditions,  but  the  study 
did  ascertain  that  the  fixed  joint  configuration  is  most  likely  the  best  overall  choice 
for  a  joined- wing  aircraft. 

In  another  journal  article  [3],  Gallman  and  Kroo,  evaluated  the  weight  opti¬ 
mization  method  and  the  fully  stressed  design  method  to  compare  the  direct  operat¬ 
ing  cost  of  conventional  wing  with  joined- wing  design.  They  determined  that  using 
a  fully  stressed  design  (FSD)  based  on  non-linear  analysis  produces  adequate  results 
with  much  less  computational  time  than  a  weight  optimization  method.  They  also 
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observed  that  adding  a  buckling  criterion  to  the  design  caused  the  joined-wing  to 
have  a  higher  direct  operating  cost  than  the  conventional  aircraft.  A  FSD  was  used 
for  the  work  [2]  preceding  this  research,  and  FSD  with  non-linear  analysis  will  be 
used  throughout  the  remainder  of  this  project.  [18] 

Livne  surveyed  these  and  over  seventy  other  studies  up  through  the  year 
2000  [7]  in  a  study  of  the  insights  and  limitations  learned  from  aeroelastic  and 
aeroservoelastic  analysis  and  design  of  joined-wing  configurations.  In  Livne’s  words, 
“there  is  still  a  basic  lack  of  understanding  of  the  structural  and  aeroelastic  behavior 
and  associated  design  issues  and  failure  modes  of  such  configurations.”  [7]  He  points 
out  the  need  for  research  in  the  areas  covered  in  the  current  work: 

Additionally,  possible  interactions  on  surveillance  JW  (joined-wing) 
aircraft  between  antennas  built  into  the  composite  lifting  surfaces  and 
the  structure  itself  should  be  studied.  Such  interactions  involve  both  the 
effects  on  stiffness,  damping  and  inertia  of  the  integrated,  multipurpose 
structure,  and  well  as  effects  of  antenna  power  consumption  and  heating. 

Strict  stiffness  requirements,  imposed  to  meet  precise  antenna  shapes  and 
required  electromagnetic  performance  can  also  add  to  the  complexity  of 
the  MDO  (multi-disciplinary  optimization)  problem.  [7] 

Livne  also  pointed  out  that  the  maneuver  load  factor  and  type  of  joint  seem  to  be 
the  most  important  factors  of  aeroelastic  behavior,  although  the  structural  design  of 
the  lifting  surfaces  was  also  important. 

M.  Blair  and  R.  Canfield  used  a  FSD  model  of  a  joined  wing  in  their  work  [2], 
which  laid  the  foundation  for  the  current  study.  Their  efforts  focused  on  creating 
a  design  process  that  generates  reliable  weight  estimations  of  joined-wing  aircraft. 
They  wanted  a  design  that  simultaneously  satisfies  the  range  requirements,  equi¬ 
librium  in  lift  and  pitching  moment,  and  stress  in  static  aeroelastic  equilibrium. 
The  structural  and  aerodynamic  interface  and  non-linear  trim  were  automated  us¬ 
ing  the  Adaptive  Modeling  Language  (AML)  to  link  ASTROS,  NASTRAN  and 
PanAir  [22,  12,  8,  16]  into  a  common  design  environment  called  the  Air  Vehicles 

Technology  Integration  Environment  (AVTIE)  [2],  A  designer  drives  the  analysis 
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with  an  AML  interface.  PanAir  is  used  for  the  aerodynamic  analysis.  It  calculates 
the  aerodynamic  load  on  the  wing  using  a  panel  method.  MSC.NASTRAN  is  used  to 
investigate  buckling  in  the  wing  and  analyze  the  geometric  nonlinearity  of  the  model. 
The  authors  concluded  that  geometric  nonlinearity  is  an  important  design  consider¬ 
ation  for  a  joined-wing  aircraft  and  should  be  included  in  future  models.  An  effort  is 
currently  underway  to  incorporate  this  nonlinearity  into  the  design  procedure.  [18] 

2.3  Integrating  Antennas  into  the  Joined-Wing  Structure 

Very  little  work  has  been  documented  with  respect  to  integrating  the  antennas 
into  the  structure  of  a  joined-wing  aircraft.  One  relevant  study  investigated  the  use 
of  a  technology  called  Active  Aeroelastic  Wing  in  a  joined-wing  design  to  minimize 
the  deformations  of  embedded  antennas.  [17]  In  this  study,  Reich,  Raveh,  and  Zink 
concurrently  used  six  control  surfaces,  as  shown  in  Figure  2.2,  to  minimize  wing 
deformation  and  trim  the  aircraft  for  1-g  level  flight.  While  this  is  a  novel  approach 
to  minimize  the  deformation  of  structurally  integrated  antennas,  it  is  important  to 
note  that  the  antennas  were  strictly  modeled  as  non-structural  masses  placed  in  the 
front  and  aft  wings.  The  antennas  were  not  load-bearing  members  of  the  structure, 
and  therefore,  the  stiffness  properties  of  the  antennas  were  not  incorporated  into  the 
design. 

The  results  of  this  study  show  that  the  Active  Aeroelastic  Wing  technology 
is  a  feasible  concept  for  minimizing  the  wing  deflection  and  thereby  minimizing 
the  antenna  deformation.  This  will  help  to  keep  the  degradation  of  the  antenna 
performance  to  a  minimum.  Their  research  also  points  out  that  the  most  outboard 
and  inboard,  trailing  edge  control  surfaces  on  the  front  wing  appear  to  be  the  most 
effective  for  trimming  the  aircraft  and  “all  three  leading  edge  surfaces  seem  to  be 
impractical  trim  devices.” 
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Figure  2.2  Description  of  control  surfaces  utilized  for  Active  Aeroelastic  Wing  [17] 
2-4  Conformal  Antenna  Theory 

Virtually  every  text  written  on  conformal  arrays  follows  the  same  theoretical 
progression.  Linear  arrays  are  first  established  and  characterized.  Then  groups  of 
linear  arrays  are  formed  to  make  a  planar  array,  and  finally,  a  linear  or  planar  array  is 
wrapped  around  the  surface  of  a  cylinder  to  make  a  conformal  array.  W.  H.  Kummer 
follows  this  same  logic  in  his  paper  [5]  on  basic  array  theory.  He  defines  the  array 
variations  for  his  paper  as  follows: 


Linear  arrays  consist  of  equally  spaced  elemental  radiators  laid  out  in 
a  straight  line,  while  two-dimensional  planar  arrays  consist  of  radiators 
oriented  on  a  geometric  grid  in  a  plane.  Rectangular  arrays  may  be 
thought  as  a  set  of  linear  arrays  placed  next  to  each  other,  equally  spaced, 
forming  the  two-dimensional  array.  A  linear  array  may  also  be  wrapped 
around  a  curved  surface,  usually  a  circle  or  a  cylinder.  Two-dimensional 
arrays  can  be  formed  by  replicating  these  linear  arrays  along  generatrices 
of  cylinders  and  cones,  or  wrapping  them  on  spheres.  [5] 
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Kummer  makes  common  simplifying  assumptions  in  his  analysis  including  a 
large  distance  from  the  source  to  the  far-held  region,  thereby  allowing  a  parallel  ray 
approximation  of  the  vectors  pointing  from  each  element  to  the  far-held  observation 
point.  He  also  assumes  that  the  elements  are  discrete  with  no  coupling  between 
elements,  and  they  are  arrayed  in  a  repetitive  grid  with  a  constant  distance  between 
elements.  Thus,  he  defines  the  electric  held  of  the  nth  element  in  any  array  as: 

En(x,  y,  z )  =  (e~3kr/ r)Aen(x',  y\  z')In  exp {jkr'n  cos£n).  (2.1) 

The  r  term  is  the  distance  from  the  source  to  the  observation  point  in  the  far-held 
region,  and  k  is  the  wave  propagation  constant,  which  is  2tt/\  in  free-space  with  A 
being  the  wavelength.  Ae  is  a  term  that  describes  the  electric  held  of  one  element 
by  its  polarization  and  orientation,  In  is  the  current  applied  to  the  nth  element,  and 
the  phase  excitation  is  described  by  the  remaining  exponential  term.  This  equation 
is  summed  over  all  of  the  array  elements,  and  the  final  equation  describing  any  linear 
array  consisting  of  n  elements  is  given  as: 

E(x,  y,  z)  =  (. e~3kr/r )  ^  AeJn  exp (jkr'n  cos£n).  (2.2) 

The  portion  of  this  expression  remaining  in  the  summation  is  referred  to  as 
the  array  factor,  Aa.  It  is  based  on  the  geometry  of  the  array  and  not  the  elements 
themselves.  Therefore  different  elements  arrayed  in  the  same  way  can  be  described 
using  the  same  array  factor. 

Kummer  then  gives  a  discussion  of  linear  array  pattern  analysis  and  synthesis 

by  examining  a  four-element  array  aligned  along  the  z-axis.  One  point  of  interest  in 

this  section  is  how  the  beam  of  the  array  is  scanned  by  adding  a  phase  progression,  a 

,  to  the  element  current  making  In  =  lQne~^a^n~l\  “At  the  peak  of  the  main  beam, 

9  =  90 ,  and  the  exponent  equals  zero,  thus  a  =  kd  cos  90,”  [5]  where  9  is  the  elevation 

angle  of  the  main  beam.  The  array  factor  then  becomes 
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(2.3) 


N 

Aa  =  ^2  ton  exp  (j  ( n  -  1)  [ kd  (cos  9  -  cos  0O)]). 

n= 1 

In  this  equation  the  value  d  is  the  uniform  spacing  between  elements. 

Linear  arrays  are  then  expanded  into  two-dimensional  arrays  by  aligning  linear 
arrays  next  to  one  another.  A  second  summation  of  m  elements  is  then  required  in 
the  array  factor  for  a  planar  array.  The  same  logic  for  beam  scanning  as  described 
previously  is  again  applied  to  the  planar  array  to  allow  scanning  in  this  second 
dimension.  Kummer  pointed  out  that  there  is  an  angular  limit  to  the  amount  that 
a  planar  array  can  effectively  be  scanned  without  significant  losses.  This  limit  is 
typically  within  70°  of  broadside  (the  vector  normal  to  the  plane).  [5] 

Conformal  arrays  are  formed  by  wrapping  linear  or  planar  arrays  around  the 
surface  of  shapes  such  as  cylinders,  cones,  and  spheres.  One  of  the  main  difficulties 
in  analyzing  conformal  arrays  is  that  the  beam  patterns  of  the  elements  can  point  in 
various  directions.  This  means  that  the  term  describing  the  element  pattern  can  no 
longer  be  factored  out  of  the  summations.  An  array  of  elements  located  on  an  arc  of 
a  circle,  as  shown  in  Figure  2.3,  is  examined.  This  array  is  projected  onto  a  line  that 
is  perpendicular  to  the  pointing  direction  of  the  main  beam.  This  equivalent  array 
is  symmetric  about  the  geometric  center  of  the  elements,  but  the  spacing  between 
each  element  is  no  longer  constant. 

Rummer’s  analysis  of  this  configuration  shows  the  development  of  large  grating  lobes 
at  large  scan  angles.  He  concluded  that  linear/planar  array  techniques  are  not  ade¬ 
quate  for  synthesizing  conformal  arrays. 

Balanis  develops  his  textbook  [1]  theory  on  arrays  similar  to  that  of  Kummer. 
Balanis  further  elaborates  on  the  phasing  required  for  an  end-fire  configuration, 
where  virtually  all  of  the  radiation  is  directed  along  the  axis  of  the  array.  He  chose 
to  use  /3  instead  of  a  to  represent  the  interelement  phase  progression,  but  otherwise 
the  equations  are  the  same.  The  total  phase  of  a  linear  array,  ^  =  kd  cos(0)  +  /3, 
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Figure  2.3  Array  conforming  to  a  circular  arc  projected  onto  a  line  forming  an 
equivalent  linear  array  with  unequal  spacing 

where  f3  is  obtained  by  setting  -0  =  0  with  9  then  becoming  90.  It  is  important  to 
note  that  this  development  allows  (3  to  be  either  positive  or  negative,  depending  on 
the  desired  angular  location  of  the  maximum  amplitude,  so  it  is  important  to  use  the 
correct  value  of  d0  when  steering  the  beam.  Therefore,  for  an  end-fire  configuration, 
90  is  set  equal  to  either  0°  or  180°.  One  can  also  see  that  by  progressively  changing 
/3,  the  direction  of  the  maximum  can  be  scanned  between  0°  and  180°. 

Another  parameter  that  Balanis  explains  has  a  large  effect  on  the  radiation 
pattern  is  the  spacing  between  elements,  d. 


If  the  element  separation  is  d  =  A/2,  end-fire  radiation  exists  in  both 
directions  (9  =  0°  and  9  =  180°).  If  the  element  spacing  is  a  multiple  of  a 
wavelength  (d  =  nX,  n  =  1,  2,  3,  . . . ),  then  in  addition  to  having  end-fire 
radiation  in  both  directions,  there  also  exist  maxima  in  the  broadside 
directions.  Thus  for  d  =  nX,  where  n  =  1,  2,  3,  . . . ,  there  exist  four 
maxima;  two  in  the  broadside  directions  and  two  along  the  axis  of  the 
array.  To  have  only  one  end-fire  maximum  and  to  avoid  any  grating 
lobes,  the  maximum  spacing  between  the  elements  should  be  less  than 
dmax  <  A/2.  [1] 
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This  methodology  is  easily  extended  for  a  rectangular  array.  In  the  XY-plane,  the 
phase  is  separated  into  two  components,  0 x  and  0^,  where  each  takes  into  account 
the  spacing  and  progressive  phasing  in  their  respective  directions,  given  by 

0a,  =  kdx  sin  9  cos  0  +  j3x  (2.4a) 

'ijjy  =  kdy  sin  d  sin  0  +  /3y.  (2.4b) 


Then  for  an  m  x  n  array  with  a  uniform  current,  Jo,  the  array  factor  (AF)  is  com¬ 
prised  of  two  distinct  summations  as: 


M 

AF  =  Jc  ■  £  gj{m—l)(kdx  sin  0  cos  (j)-\-(3x) 

m—  1  n= 1 


N 

^  ^  ^j(n—l)(kdy  sin  6  sin  </>+f3y) 


(2.5) 


Rewriting  this  equation  in  terms  of  sine  functions  and  normalizing  by  the  constant 
current  magnitude  the  equation  is  then  given  as 


AFn(d,0) 


1  sin  (f  0a,) 
M  sin  (^f) 


1  sin  (f  ij}y)  \ 

N  si0f)  s 


(2.6) 


Even  though  this  equation  was  developed  for  isotropic  point  sources,  it  is  useful 
because  the  total  radiation  field  of  an  array  of  any  type  of  elements  can  be  found 
by  multiplying  the  array  factor  of  the  isotropic  source  times  the  radiation  pattern  of 
the  actual  source.  This  pattern  multiplication  can  be  very  useful  for  geometrically 
simple  arrays,  but  every  element  in  the  array  must  be  identical. 

In  the  Antenna  Engineering  Handbook  [9],  R.  J.  Mailloux  makes  one  other 
observation  about  conformal  arrays  that  could  be  useful  for  certain  applications  of 
structurally  embedded  antennas.  He  points  out  that  if  the  length  of  the  array  is  small 
relative  to  the  radius  of  curvature  of  the  surface  to  which  the  array  is  conforming, 
then  the  array  may  be  considered  flat  in  order  to  simplify  the  problem.  This  means 
that  the  array  may  be  treated  as  if  it  was  planar,  and  the  methods  described  above 
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can  be  utilized.  He  does  not  further  quantify  what  ratio  of  length  to  radius  is 
appropriate  for  this  to  be  a  reasonable  approximation.  The  type  of  airfoil  used  in 
the  aircraft  design  would  determine  the  practicality  of  this  assumption. 
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III.  Methodology  and  Modeling 


3.1  Antenna  Structural  Modeling 

A  simple  structural  model  of  the  sensors  was  created,  as  shown  in  Figure  1.3, 
to  incorporate  the  different  layers  of  materials  that  make  up  the  antenna’s  structure. 
The  model  has  five  layers  that  represent  the  main  load  bearing  portions  of  the  an¬ 
tenna.  These  finite  element  models  were  incorporated  into  the  wing  structure  where 
the  arrays  are  to  be  located,  replacing  the  aluminum  that  was  used  for  initial  opti¬ 
mization  efforts  [18] .  This  allows  a  more  accurate  model  of  the  wing  structure  to  be 
generated.  This  also  affects  the  overall  stiffness  of  the  structure,  which  in  turn  affects 
the  amount  of  deflection  in  the  wing.  The  simplified  model  consists  of  an  electro- 
magnetically  transparent  material  called  Astroquartz,  honeycomb  core  material,  and 
graphite  epoxy.  The  top,  outermost  layer  of  the  element  is  Astoquartz.  This  mate¬ 
rial  covers  the  radiating  parts  of  the  antenna,  protecting  them  from  the  atmosphere, 
while  also  allowing  most  of  the  antenna’s  energy  to  pass  through  with  minimal  re¬ 
sistance.  There  are  two  layers  of  graphite  epoxy  that  act  as  the  main  load-bearing 
components  of  the  antenna  element.  Both  the  Astroquartz  and  graphite  epoxy  layers 
were  constructed  of  0°,  ±45°,  and  90°  plies.  Finally,  two  layers  of  honeycomb  core 
material  act  as  spacers  between  the  other  materials.  It  is  in  these  two  layers  that 
the  launchers,  combiners,  and  other  antenna  components  are  located.  Since  these 
elements  do  not  carry  any  of  the  loads,  they  can  be  ignored  in  the  structural  analysis. 
The  thickness  of  each  layer  was  established  by  Northrop  Grumman  [15]  based  on  the 
required  separation  between  various  antenna  components  and  optimization  of  the 
structure’s  weight  and  rigidity.  The  strength  properties  of  each  layer  are  outlined  in 
Table  3.1. 

Some  additional  values  were  also  used  for  the  honeycomb  core  layers  [4],  Two 
shear  moduli  were  given  as  15.0  ksi  and  8.0  ksi,  corresponding  to  the  L  and  W 
directions,  respectively,  as  shown  in  Figure  3.1.  Also,  the  density  of  the  chosen 
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Table  3.1  Mechanical  properties  of  Astroquartz,  graphite  epoxy,  and  honeycomb 


core  [15,  4] 


Material 

Ex  [ Msi 

Ey  [Msi 

^xy 

Gxy  [ Msi 

tply 

Astroquartz 

6.80 

1.34 

0.36 

0.72 

0.0055 

Graphite  Epoxy 

22.13 

2.15 

0.3 

0.6 

0.0056 

Honeycomb  Core 

0.07 

0.07 

— 

0.0004 

0.9/0. 5 

material  was  4.5  lb/ in3,  and  the  compressive  strength  was  560  ksi,  which  was  assumed 
to  be  the  same  in  both  the  x  and  y  directions. 


w 


A 


L 


Figure  3.1  L  and  W  directions  of  honeycomb  core  material 

To  test  the  validity  of  this  model,  buckling  analysis  was  performed  using  NAS- 
TRAN  on  a  rectangular  mesh  of  these  elements,  which  represented  a  typical  wing 
skin  panel.  The  panel  was  33  inches  in  the  span  direction  and  44  inches  in  the  chord 
direction.  The  initial  mesh  contained  22  elements  in  each  direction.  A  compressive 
load  of  12,000  lbs. /in  was  applied  to  the  longer  sides  of  the  panel.  This  represented  a 
spanwise  load  carried  along  the  wing.  The  panel  was  clamped  along  the  two  longer, 
chordwise  edges  and  simply  supported  along  the  other  two  spanwise  edges.  The  first 
buckling  mode  eigenvalue  and  shape  were  obtained,  and  the  mesh  was  then  refined 
to  test  for  convergence.  The  mode  shapes  were  plotted  using  I-DEAS  [20].  These  re¬ 
sults  were  compared  to  those  generated  by  Northrop  Grumman  [15],  and  the  results 
are  shown  in  Chapter  4,  Section  1. 
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3.2  Antenna  Electromagnetic  Model 

The  conformal  load  bearing  antenna  element  is  a  complicated  structure  that 
would  have  been  extremely  difficult  to  model  as  a  single  entity.  Therefore,  it  was 
decided  that  simple  dipole  (wire)  antenna  elements  would  be  used  to  clectromag- 
netically  approximate  the  actual  elements.  Dipole  antennas  are  one  of  the  simplest 
antenna  structures,  and  are  widely  used  in  the  development  of  linear  and  planar 
array  theory.  All  of  the  dipoles  were  assumed  to  be  identical,  perfect  conductors 
with  an  infinitesimally  small  diameter. 

The  CLAS  elements  were  designed  to  produce  a  horizontally  polarized  electric 
field.  To  represent  this  polarization,  the  dipoles  were  oriented  parallel  to  the  wing 
surface  with  5  elements  arrayed  in  the  chordwise  direction  of  the  wing,  conforming 
to  the  shape  of  the  airfoil,  and  40  elements  in  the  spanwise  direction.  The  0  compo¬ 
nent  of  the  electromagnetic  field  (see  Equation  3.21)  then  describes  the  horizontally 
polarized  portion  of  the  total  field. 

A  baseline  configuration  of  the  array  was  generated  to  correspond  to  the  an¬ 
tenna  orientation  when  the  wing  is  not  deformed.  In  this  case  the  wing  is  assumed 
to  be  straight  and  at  a  fixed  dihedral  angle.  The  antenna  array  is  positioned  inboard 
of  the  wing  joint  and  as  close  to  the  joint  as  possible  considering  the  current  30° 
wing  sweep.  This  placement  is  shown  in  Figure  3.2. 

The  mesh  used  to  generate  the  coordinate  data  hie  was  created  by  intersecting 
the  wing  surface  with  2N  +  1  planes  in  the  spanwise  direction,  and  2 M  +  1  planes 
in  the  chordwise  direction,  where  N  and  M  describe  the  number  of  elements  in  the 
spanwise  and  chordwise  directions,  respectively.  The  factor  of  two  was  used  so  that 
the  coordinates  of  the  center  point  of  each  element  were  obtained,  as  well  as  the 
coordinates  for  the  point  that  bisected  each  side.  Therefore,  nine  nodes  described 
the  location  and  orientation  of  each  element.  A  section  of  a  generic  planar  antenna 
with  square  elements  is  shown  in  Figure  3.3  to  illustrate  the  nine  node  configuration. 
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Figure  3.2 


Antenna  position  on  front  section  of  joined  wing 


The  two  points  in  the  middle  of  each  element’s  chordwise  side  determined  the 
length  of  the  dipoles  utilized  in  this  analysis.  For  an  antenna  frequency  of  430  MHz, 
this  length  was  approximately  half  of  a  wavelength.  The  center  points  of  the  actual 
elements  correspond  to  the  center  points  of  the  dipoles.  Initially,  five  dipoles  were 
used  to  represent  the  five  CLAS  elements  in  the  chordwise  direction,  but  the  final 
configuration  of  the  dipoles  utilized  10  dipoles  to  more  closely  model  the  effective 
spacing  of  the  CLAS  elements.  The  reasoning  behind  this  change  is  further  explained 
in  the  next  chapter. 

Two  other  versions  of  the  array  were  used  besides  the  10x40  front  wing  array. 
First  the  array  size  was  increased  to  10x60  elements  on  the  front  wing.  This  version 
is  shown  in  Figure  3.4.  The  second  version  placed  a  10x60  element  array  in  the  aft 
wing  section,  which  is  shown  in  Figure  3.5. 

End-fire  phasing  was  applied  to  steer  the  main  beam  straight  out  from  the  lead¬ 
ing  edge  and  in  the  plane  of  the  wing.  Once  the  baseline  pattern  was  established  for 
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Center  Points  of  Dipoles 


Dipole 


Figure  3.3  Nodes  on  square,  S-CLAS  element  representing  where  coordinates  are 
required  for  analytical  development  of  the  dipole  array 


Figure  3.4  10x60  element  array  located  on  the  front  wing  section 

the  undeformed  case,  specific  load  configurations  were  applied  to  the  lifting  surfaces. 
Wing  deformations  were  then  generated  on  a  fully  stressed  design  [18]  by  an  inte¬ 
grated  software  environment  called  AVT1E  [2],  which  utilizes  the  Adaptive  Modeling 
Language  (AML),  MSC.NASTRAN  and  PanAir  [22,  8,  16].  These  deformations  were 
used  to  determine  the  new  location  and  slope  of  each  element  of  the  array.  Figure  1.4 
shows  the  typical  deformed  shape  of  the  array  compared  to  the  undeformed  shape. 
A  new  beam  pattern  was  then  generated.  This  process  was  repeated  for  various 
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Figure  3.5  10x60  element  array  located  on  the  aft  wing  section 

load  configurations.  The  results  were  compared  to  the  undeformed  wing  results  to 
determine  the  beam  pointing  errors  produced  by  the  wing  deformations. 

The  radii  of  curvature  of  the  top  and  bottom  surfaces  of  the  wing  were  also 
checked  to  see  whether  Mailloux’s  flat  array  approximation  could  be  applied.  The 
local  radius  of  curvature  of  a  point  on  the  airfoil  was  determined  based  on  its  location 
relative  to  two  adjacent  points.  These  three  points  defined  an  arc  of  a  circle  and  the 
radius  of  that  circle  was  then  calculated  by  symbolically  solving  a  system  of  three 
simultaneous  equations.  The  three  equations  were 


(aq  -  x0 )2  +  (yi  -  y0)2  -  r2 

=  0 

(3.1a) 

(x2  -  x0 )2  +  (y2  -  y0)2  -  r2 

=  0 

(3.1b) 

{x3  -  xo)2  +  {V3  ~  Vo)2  -  r2 

=  0 

(3.1c) 

where  Xq  and  y0  were  the  coordinates  of  the  center  of  the  circle  with  radius,  r,  formed 
by  the  three  points  with  coordinates,  aq  and  y*,  i  —  1,2,  3.  The  LRN-1015  airfoil  is 
shown  in  Figure  3.6. 
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p - 2.5  m - p 

Figure  3.6  LRN-1015  airfoil 

A  frequency  of  430  MHz  corresponds  to  a  wavelength  of  27.5  inches.  Therefore, 
the  chordwise  length  of  the  array  was  approximately  2.5  wavelengths.  A  radius 
was  assumed  to  be  much  larger  than  this  array  length  if  it  was  over  10  times  the 
wavelength. 

3.3  Array  Theory 

Developing  the  theory  for  this  array  began  by  defining  the  geometry  and  ori¬ 
entation  of  the  array.  The  implementation  was  general  enough  to  analyze  any  array. 
This  allowed  the  same  code  to  be  utilized  for  the  array  under  any  load  configuration. 
The  coordinates  were  given  in  terms  of  a  global  wing  reference  frame,  where  the 
origin  was  located  at  the  leading-edge  of  the  root  of  the  wing.  The  orientation  of 
the  x  and  y  axes  with  respect  to  the  wing  are  shown  in  Figure  3.7. 


Figure  3.7  Global  wing  coordinate  reference  frame 
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The  reference  frame  utilized  in  the  array  development  was  a  clockwise  rotation 
of  the  x  and  y  axes  in  the  global  reference  frame  corresponding  to  the  wing  sweep. 
The  origin  was  also  translated  to  coincide  with  the  center  point  of  the  dipole  situated 
in  the  inboard,  leading-edge  corner  of  the  array,  which  was  defined  to  be  the  original 
element  that  was  arrayed  to  form  the  entire  sensor  suite.  This  reference  frame  is 
shown  in  Figure  3.8.  Two  chordwise  arrays  of  dipoles  were  drawn  into  this  figure 
to  illustrate  that  the  origin  of  the  array  reference  frame  coincided  with  the  center 
point  of  the  inboard,  leading-edge  dipole,  and  the  y- axis  lined  up  parallel  to  the 
same  inboard,  leading-edge  dipole.  Therefore,  wherever  that  particular  dipole  was 
situated  on  the  wing  defined  the  location  and  orientation  of  the  array  coordinate 
system. 


Figure  3.8  Antenna  coordinate  reference  frame  with  two  inboard  chordwise  arrays 
drawn  in  for  reference 

A  vector,  fo,  was  formed  to  point  from  the  origin  to  any  point  on  the  dipole 
centered  at  the  origin  of  the  coordinate  reference  frame,  given  by 


fo  =  yy 


(3.2) 
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The  value  y  corresponds  to  a  position  along  the  length  of  each  dipole.  For  the 
undeformed  wing  case,  all  of  the  dipoles  are  assumed  to  be  aligned  parallel  to  the  y 
direction. 

Next,  a  vector  was  formed  that  points  from  the  center  of  the  original  element 
to  the  center  point  of  the  A;th  element  in  one  chordwise  array, 

rk  =  (■ xk  ~  x0)x  +  (yk  -  y0)y  +  (zk  -  z0)z.  (3.3) 

Since  the  original  dipole  was  assumed  to  be  centered  at  the  origin  of  the  antenna 
reference  system,  this  equation  can  be  simplified  to 


rk  =  xkx  +  yky  +  zkz.  (3.4) 

Equations  3.2  and  3.4  were  then  summed  to  give  a  new  vector  that  points  from  the 
origin  to  any  point  on  the  kt\i  dipole, 

ffe  =  r0  +  rk.  (3.5) 

A  third  vector  was  also  defined.  This  was  a  unit  vector  that  pointed  in  the 
direction  of  an  observation  point  in  the  far-held  of  the  antenna.  This  vector  is  defined 
as  follows: 


f(9,  0)  —  x  sin  6  cos  0  +  y  sin  6  sin  0  +  z  cos  9  (3.6) 

where  elevation,  9,  and  azimuth,  0,  are  defined  according  to  Figure  3.9. 

For  evaluation  purposes,  the  position  of  the  observation  point  was  varied  to 

give  three  “great  circle”  cuts.  In  each  cut  one  angle  was  held  fixed  while  the  other 

was  varied  from  0  to  2tc.  To  get  an  elevation  cut,  0  was  set  equal  to  180°.  For  an 

azimuth  cut,  6  was  set  to  90°.  Finally,  for  the  third  cross-sectional  cut,  0  equaled 
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z 


Figure  3.9  Array  coordinate  system 

90°  and  6  was  varied.  These  three  cuts  are  the  most  common  for  defining  radiation 
patterns. 

The  phase  of  the  array  in  free  space  was  then  described  by  the  expression  e3^'r  1 
where  fi  is  again  the  wave  propagation  constant  for  lossless  space,  f3=u^/~jle.  In  this 
equation  for  (3,  n  and  e  are  the  permeability  and  permittivity,  respectively,  of  the 
material  surrounding  the  antenna,  and  the  angular  frequency,  ui,  is  defined  as  tv  = 
27 r/,  where  /  is  the  operating  frequency  of  the  antenna.  Substituting  Equation  3.5 
into  the  dot  product,  the  phase  then  becomes  e-?/3(r°'r+rfc'rb 

A  cosine  current  distribution  was  applied  to  each  dipole.  This  distribution 
is  defined  by  a  maximum  value  of  1  at  the  center  of  the  dipole,  going  to  zero  at 
each  end.  The  cosine  function  was  chosen  because  it  had  no  discontinuities  along 
the  length  of  the  dipole,  and  it  matched  the  boundary  conditions  of  the  current 
going  to  zero  at  both  ends  of  the  dipole.  The  current  was  a  complex  vector  quantity 
(represented  by  a  bar  under  the  symbol)  defined  below: 

Jr(0,  (p,uj)  =  ^2  Vak^5k  cos  (jrVk)  (3.7) 

k 

where  cq.  was  a  weighting  coefficient  that  could  be  used  to  scale  the  magnitude  of 
the  current  on  the  kth  element,  and  L  is  the  length  of  the  dipoles,  which  was 
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The  free  space  wavelength  A  is  defined  as  A=y.  The  coefficient  was  simply  set 
to  a  value  of  one  for  all  elements  for  the  purpose  of  this  study,  but  it  was  important 
to  incorporate  into  the  derivation  for  future  work  to  optimize  the  radiation  pattern. 
The  interelement  phasing,  5k,  was  applied  to  each  element  to  steer  the  main  beam 
in  the  desired  direction.  The  subscript  R  denotes  that  this  equation  was  derived  for 
a  rectangular  coordinate  system,  and  it  is  a  function  of  elevation,  9 ,  azimuth,  0,  and 
the  angular  frequency,  c o,  which  is  2n  times  the  radial  frequency  of  the  system. 

This  equation  was  then  expressed  in  term  of  the  frequency  domain  through  the 
use  of  a  Fourier  Transform.  The  equation  of  the  transform  was  as  follows: 


Jr(9,(J),uj) 


J  R^9(ro-f+rk-r) 


(3.8) 


An  assumption  was  made  here  that  loss  is  ignored  in  the  integral.  This  assumption 
is  valid  because  the  loss  has  no  effect  on  frequency,  and  it  is  transverse  to  the  main 
direction  of  propagation  and  loss,  which  is  in  the  direction  of  r. 

Taking  the  derivative  of  Equation  3.5  and  assuming  that  the  kt h  element  is  at 
a  fixed  point  in  space,  then 


dr*k  =  dr0.  (3.9) 

Noting  also  that  the  derivative  of  r0  with  respect  to  yk  is  simply  dyk,  a  change  of 
variables  was  performed  on  Equation  3.8,  and  the  equation  for  J_R  became 

JR{e,M  =  J  lRe^+^dyk.  (3.10) 

Vk 

The  constant  terms  were  then  removed  from  the  integral,  and  the  final  form  of  this 
equation  before  the  transform  is  performed  is  shown  below: 
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=  v 


y^akeJ"Ke 


,jSk0j/3r-rk 


2  7T 

cos  I  —yk  I  e 


jP^°dyk. 


(3.11) 


Vk 


A  detailed  explanation  of  the  Fourier  Transform  is  given  in  Appendix  B.  A 
closed  form  approximation  was  used  for  this  study  to  simplify  the  transform.  There¬ 
fore,  the  integral  portion  of  Equation  3.11  results  in  the  following  expression: 


J  cos  eJpr'r°dyk 

Vk 


|  cos(w) 


(3.12) 


where  L  is  the  length  of  the  dipole,  and  u  for  a  ^-directed  current  is  defined  as 


u  = 


—  sin  9  sin  6. 
2 


(3.13) 


Substituting  Equations  3.12  and  3.13  back  into  Equation  3.11,  the  transformed  cur¬ 
rent  became 


where 


=  v 


V  n,  pjSkJ0r-rk 

£  cos  |  sin(0)  sin(0) 

k 

.1  -  (£sin(0)  sin(0))2_ 

(3.14) 


r  ■  rk  =  ( xk  —  Xq)  sin  9  cos  4>  +  (yk  —  2/o)  sin  9  sin  0  +  (zk  —  z0)  cos  9  (3.15) 


A  vector  potential,  ARl  in  rectangular  coordinates  was  then  created.  This  is 
common  practice  when  solving  Maxwell’s  Equations,  because  it  allows  a  difficult 
integration  to  be  broken  up  into  a  simpler  integration  and  differentiation  steps  as 
shown  in  Figure  3.10.  The  vector  potential  in  a  lossless  medium  is  defined  as 
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Figure  3.10  Pictorial  representation  of  the  use  of  vector  potentials  to  find  the 
electromagnetic  fields 

Ar  =  fi0^(r,u)  JR.  (3.16) 

where  fi0  is  the  permeability  of  free-space  (47rxl0-7),  and  4/(r,  uj)  is  the  radial  com¬ 
ponent  of  Green’s  Function  in  free-space  defined  as  follows: 


T  (r,  u) 


e-JPr 
47T  r 


(3.17) 


where  r  is  the  distance  from  the  source  to  the  far- field  observation  point. 

The  vector  potential  was  then  transformed  from  rectangular  to  spherical  coor¬ 
dinates  by  use  of  the  transformation  matrix,  Tg. 


sin  9  cos  0  sin  9  sin  0  cos  9 

=R 

Tg  =  cos  9  cos  0  cos  9  sin  0  —  sin  6  (3.18) 

—  sin  0  cos  0  0 

Since  the  current  was  assumed  to  be  oriented  entirely  along  the  y-axis,  the 
x  and  z  components  of  AR  were  set  equal  to  zero,  and  Ay  was  multiplied  times 
the  transformation  matrix.  Also,  the  radial  component  of  the  vector  potential  in 
spherical  form  was  discarded  based  on  the  far-field  assumption  that  there  is  no 

variation  in  the  radial  direction.  This  resulted  in  the  following  equation: 

3-13 


(3.19) 


—S  ~  Ts  '  Ar- 


The  far-field  region  E- field  in  spherical  coordinates  minus  the  radial  component 
was  then  defined  as 


Es  =  -ju  [As  -  r Ar]  .  (3.20) 

Expanding  this,  the  equation  for  the  E-field  was  as  follows: 


e~iPr 

Es  = 


Y  akej5kej0^k 

k 


-  cos  r 4  sin  9  sin  0 1 

r  -| 

it  L  2  “  J 

,  Y  <  v  2 

6  cos  9  sin  0  +  0  cos  0 

1  —  ( ^  sin  9  sin  0) 

(3.21) 


where,  for  horizontal  dipoles,  the  9  component  corresponds  to  the  vertically  polarized 
portion  of  the  equation,  and  the  0  component  corresponds  to  horizontal  polarization. 
Since  this  study  was  interested  in  emulating  a  horizontally  polarized  field,  only 
was  used.  Therefore  the  plots  that  were  generated  correspond  to: 


-jcj/do 


e~jf3r 

4-7T  r 


Y  akejSkejpf  Pk 

.  k 


^  cos  sin  9  sin  0 
1  —  (^  sin  9  sin  0) 2 


COS  0. 


(3.22) 


The  MATLAB  code  used  to  plot  the  three  Great  Circle  cuts  of  Equation  3.22 
breaks  the  sum  on  k  into  two  summations:  one  to  sum  the  elements  in  the  chordwise 
direction,  forming  a  chordwise  array,  and  then  a  second  to  sum  up  all  of  the  chordwise 
arrays.  This  summing  process  was  a  simplification  that  originates  from  the  far-field 
planar  array  theory. 

Thus  far,  the  equations  were  developed  assuming  that  the  wing  was  unde¬ 
formed.  Therefore  it  was  assumed  that  each  new  dipole  orientation  and  location 
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was  strictly  due  to  a  translation  of  the  element  from  the  origin.  For  a  deformed 
wing  this  assumption  was  no  longer  valid,  and  the  amount  of  rotation  experienced 
by  each  element  had  to  be  calculated  and  incorporated  into  the  radiation  pattern 
equations.  To  accomplish  this,  the  same  data  hie  described  in  Section  3.2,  page  3-3, 
was  used  in  conjunction  with  a  central  difference  formula  to  calculate  the  slope  of 
each  dipole.  The  central  difference  formulation  utilizes  the  y  and  z  coordinates  of 
the  data  points.  A  forward  difference  method  was  used  to  characterize  the  slope  of 
the  elements  along  the  inboard  edge  of  the  array,  and  a  backward  difference  was  used 
along  the  outboard  edge.  These  equations  are  defined  as  follows: 


slope(n,  1) 
slope  (n,  m) 


z(;n,  2)  —  z(n,  1) 
h 

z(n,  m)  —  z(n,  m 
h 


1) 


(3.23) 

(3.24) 


where  h  was  the  difference  in  y  of  adjacent  nodes,  which  equals  y(n,  2)  —  y(n,  1) 
in  Equation  3.23  and  y(n,m)  —  y(n,  m  —  1)  in  Equation  3.24.  The  indices  n  and 
m  describe  the  number  of  data  points  in  the  chordwise  and  spanwise  directions, 
respectively. 

For  the  interior  points  of  the  array,  the  following  equation  was  used: 


slope(n,  m  +  1) 


(yg SA(yH+J^y(^+2)) 

h(n,  m)  +  h(n,  m  +  1) 

1  1 


h(n,  m)  h(n,  m  +  1) 


z(n,  m  +  1). 


(3.25) 


This  new  orientation  of  each  dipole  was  then  incorporated  into  Equation  3.13.  This 
equation,  which  described  the  orientation  of  the  current  on  each  dipole,  then  became 

u(n,m)  =  —  (cos(slope(n,  m))  sin 0 sin 0  +  sm(slope(n,  m))  cos#),  (3.26) 
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and  was  used  as  before  in  the  equation  describing  the  E-field,  Equation  3.21. 

With  the  dipoles  situated  parallel  to  the  surface  of  the  wing,  any  twist  in  the 
wing  resulted  in  only  a  translation  of  the  dipole.  These  translations  were  accounted 
for  in  the  vector  r*,,  since  it  pointed  to  the  center  of  the  dipole  in  its  new  position. 
The  data  hies  used  to  generate  the  patterns  for  the  deformed  wing  cases  included 
the  trim  angle  of  attack  and  aft  wing  rigid  twist  that  correspond  to  each  load  case. 
Therefore,  the  trimmed  angle  of  attack  was  subtracted  from  the  main  beam  pointing 
angle  results  so  that  only  changes  clue  to  wing  deflection  are  present  for  the  array  in 
the  front  wing.  For  the  aft  wing  array,  the  twist  required  to  trim  for  pitch  remained 
in  the  results,  because  the  twist  is  not  constant  out  the  span  of  the  aft  wing.  For 
the  2.5g  maneuver  load  the  trim  angle  of  attack  was  0.26°,  and  the  aft  wing  twist 
was  —4.23°.  For  the  gust  load  case,  the  angle  of  attack  was  1.28°  and  the  aft  wing 
twist  was  1.08°. 

3.4  Computational  EM  Method 

For  comparison  purposes,  the  model  was  also  analyzed  with  a  numerical  elec¬ 
tromagnetic  code  that  solves  for  the  radiation  pattern  of  the  array  using  a  method 
of  moments  technique.  The  software  package  is  called  NEC-Win  Plus+™,  and  it 
is  a  scaled  down,  more  user-friendly  version  of  the  full  NEC  software  [14],  It  allows 
simple  antenna  models  to  be  built  fairly  quickly,  and  it  can  numerically  determine 
their  electromagnetic  response  for  a  given  current  excitation.  The  properties  for  the 
antenna  array  were  defined  to  have  a  frequency  of  430  MHz,  no  ground  plane,  a  char¬ 
acteristic  impedance  of  50  Ohms,  each  wire  constructed  of  11  segments,  a  current 
source  at  the  center  of  the  wire,  current  magnitude  of  1.0,  a  wire  diameter  of  0.001 
meters,  a  perfectly  conducting  wire,  and  90°  phasing  applied  from  one  element  to  the 
next  in  the  chordwise  direction.  The  electric  held  integral  equation  was  described 
in  the  previous  section,  but  there  a  closed  form  solution  was  used  to  approximate 
the  solution  of  the  Fourier  Transform  clue  to  the  simplicity  of  the  elements.  NEC- 
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Win  Plus+™  numerically  solves  these  equations  to  give  a  closer  approximation  of 
the  exact  solution.  This  process  assumes  nothing  about  the  orientation  of  one  cl¬ 
ement  with  respect  to  another.  By  using  this  software,  it  was  possible  to  see  how 
applicable  basic  array  theory  was  for  the  given  array.  It  was  preferable  to  use  the 
MATLAB  code,  because  it  was  computationally  less  intensive  and  could  be  more 
easily  incorporated  into  the  aircraft  design  environment. 

The  method  of  moments  is  a  general  technique  for  approximating  the  solution 
of  linear  operator  equations.  It  does  this  by  creating  a  system  of  simultaneous 
linear  equations  that  represent  the  linear  operator  equation.  For  example,  the  linear 
operator  equation 

Lf  =  g  (3.27) 

consists  of  the  linear  operator,  L,  the  unknown  function,  /,  and  some  known  function, 
g.  In  the  context  of  electromagnetics,  one  typically  knows  the  integral  equation 
describing  the  electric  or  magnetic  field.  This  corresponds  to  the  linear  operator, 
and  g  represents  the  known  excitation  of  the  system  (i.e.  source  currents).  The 
unknown  function,  /,  is  therefore  the  radiated  energy  due  to  the  excitation.  A  series 
expansion  is  used  to  approximate  the  unknown  function. 

f  =  YJanIn  (3.28) 

n 

In  this  equation,  an  is  a  scalar  coefficient  applied  to  the  basis  function,  In.  The  choice 
of  basis  functions  determines  the  accuracy  of  the  solution.  Therefore,  a  great  amount 
of  consideration  typically  goes  into  selecting  appropriate  basis  functions  for  a  given 
system.  The  reader  is  referred  to  the  NEC  manual  [13]  for  a  complete  derivation  of 
the  basis  functions  which  satisfy  the  boundary  conditions  for  wire  elements. 
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An  inner  product  of  two  functions  is  then  defined  as  follows: 


(x,v)  =  I  x(r)v(r)dA  (3.29) 

s 

integrating  over  the  entire  surface  of  the  system.  The  inner  product  is  taken  of  Lf 
and  a  set  of  test  functions,  wm .  The  NEC  code  chooses  these  test  functions  to  be 
delta  functions, 

wm(f)  =  8(f—  rm)  (3.30) 

where  f  points  to  the  far-held  observation  location,  and  rm  points  to  a  location  on 
the  surface  of  the  antenna.  These  points  on  the  surface  of  the  antenna  are  the  center 
points  of  each  segment  on  every  antenna.  This  point  sampling  method  was  the  main 
reason  why  this  approach  took  a  considerable  amount  of  time  and  computer  capacity. 
For  the  antenna  array  studied,  this  corresponded  to  400  elements  with  11  segments 
each,  or  4400  points  at  every  angular  step  in  a  360°  observation  range. 

The  inner  product  of  Equation  3.27  and  the  test  functions  is 

(wm,Lf)  =  (wm,g) ,  (3.31) 

where  rn  —  1, 2, 3, . . . ,  N.  Equation  3.28  is  then  substituted  into  this  equation.  The 
resulting  inner  product  is  then 

N 

^  ^  cxn  (wm)  LIn )  (wrni  g) .  (3.32) 

n=l 

This  equation  can  be  written  in  matrix  form  as 

[Z]  M  =  [G]  (3.33) 
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where  the  matricies  are  defined  as  follows: 

(w1}LIi)  (wi,  LI2)  ■■■ 

[Z\  =  (w2,LI1)  (w2,LI2)  •••  (3.34) 

Ai 

[A]  =  A2  (3.35) 

(wi  ,g) 

[G]  =  (w2,g)  ,  (3.36) 

and  the  solution  is 

[A]  =  [Z]-‘[G],  (3.37) 

In  an  electromagnetic  context,  Z  is  typically  referred  to  as  the  impedance  matrix,  G 
is  the  excitation  vector,  and  A  is  called  the  current  vector.  Based  on  this  derivation, 
Equation  3.28  then  becomes 

N 

f  =  YJAnIn  (3.38) 

n=  1 

for  the  chosen  set  of  basis  functions,  /. 

3.5  Control  Surface  Effectiveness 

An  aerodynamic  model  was  created  using  MSC.FlightLoads  [11].  The  model 
was  generated  for  a  2.5g  pull-up  maneuver  at  the  end  of  the  mission  profile  (i.e. 
landing),  therefore  sea  level  conditions  were  assumed.  This  model  consisted  of  four 
main  sections:  the  front  wing  inboard  of  the  joint,  the  aft  wing  inboard  of  the  joint, 
the  joint,  and  the  remainder  of  the  wing  outboard  of  the  joint  (wing  tip).  The 
aerodynamic  model  is  shown  in  Figure  3.11.  An  aileron  was  placed  in  the  wing  tip 
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section  of  this  model  and  used  to  trim  for  pitch  so  the  effectiveness  of  this  control 
surface  could  be  determined.  This  control  surface  was  sized  to  be  30%  of  the  chord 
length  and  half  of  the  span  length  of  the  wing  tip  section,  and  located  along  the 
trailing  edge  of  the  wing  tip  section  as  far  outboard  as  possible. 


Figure  3.11  Aerodynamic  panel  model  of  joined  wing  with  shaded  panels  showing 
location  of  aileron 

Once  the  aerodynamic  model  was  built,  FlightLoads  was  used  to  spline  it  to 
the  structural  model.  FlightLoads  then  generated  a  NASTRAN  input  hie,  and  it  was 
analyzed  to  find  the  static  aeroclastic  response  of  the  model.  From  the  NASTRAN 
output  hie,  the  coefficient  of  lift  per  unit  change  in  control  surface  deflection  (5)  of 
an  elastic  wing,  Ci,5(elastic ),  and  the  coefficient  of  lift  per  unit  change  in  5  of  a  rigid 
wing,  CLs(rigid),  were  obtained.  These  two  values  are  used  to  define  the  control 
surface  lift  effectiveness,  tjl,  as  follows: 


CLs(elastic) 
llL  CLs  (rigid) 


(3.39) 
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The  purpose  of  an  aileron  is  not  to  trim  the  aircraft  for  pitch,  so  this  effective¬ 
ness  was  expected  to  be  quite  low.  If  control  surfaces  are  not  effective  for  pitching 
the  aircraft,  though,  other  options  must  be  utilized.  This  meant  that  either  the  aft 
wing  had  to  twist  to  trim  the  aircraft  for  pitch,  or  a  control  surface  had  to  be  placed 
somewhere  inboard  of  the  joint,  which  would  limit  the  maximum  size  of  the  antenna 
arrays.  It  was  decided  that  aft  wing  twist  would  be  used.  Therefore,  the  effectiveness 
of  the  aileron  was  checked  for  its  intended  purpose,  which  was  to  roll  the  aircraft. 
NASTRAN  was  used  to  find  the  non-dimensional  roll  rate,  7^,  at  which  aileron 
reversal  occurred,  where  p  is  the  roll  rate,  b  is  the  total  span,  V  is  the  velocity, 
and  5  is  the  aileron  deflection.  The  dynamic  pressure,  q  (and  therefore  the  velocity) 
related  to  the  reversal  condition  was  also  found.  The  results  of  this  analysis  are 
shown  in  Section  4.5. 
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IV.  Results 


4.1  FEM  Results 

The  composite  antenna  elements  were  examined  by  creating  a  rectangular 
33”x44”  mesh,  similar  to  a  wing  panel,  consisting  of  these  elements.  This  was 
done  in  an  effort  to  duplicate  results  generated  in  a  similar  proprietary  effort  by 
Northrop  Grumman  [18],  thereby  validating  this  model.  Northrop  Grumman’s  three- 
dimensional  FEM  was  more  robust  than  the  two-dimensional  model  developed  for 
this  study,  but  the  model  used  gave  a  first  buckling  mode  eigenvalue  that  was  within 
acceptable  limits  of  Northrop  Grumman’s  results,  considering  the  current  fidelity  of 
the  entire  joined-wing  model. 

The  first  buckling  mode  eigenvalue  was  1.27  for  the  panel  under  a  spanwise 
compressive  load  on  the  two  longer  edges  with  all  of  the  edges  simply  supported. 
This  mode  shape  is  shown  in  Figure  4.1.  Northrop  Grumman’s  value  for  this  mode 
was  0.99.  ft  was  believed  by  Northrop  Gurmman  that  a  more  realistic  represen¬ 
tation  of  the  constraints  had  the  chordwise  edges  clamped  and  the  spanwise  edges 
simply  supported.  Therefore,  the  boundary  conditions  were  changed  to  reflect  these 
conditions.  The  eigenvalue  remained  the  same,  but  the  mode  shape  changed  to  two 
sine  waves,  as  shown  in  Figure  4.2,  instead  of  one  and  a  half  in  the  case  where 
all  of  the  edges  were  simply  supported  .  Northrop  Grumman’s  eigenvalue  for  this 
clamped/simply  supported  condition  was  1.31,  and  the  mode  shape  looked  similar 
to  Figure  4.3. 

Figure  4.3  was  achieved  for  the  model  developed  in  this  effort  by  decreasing 
the  overall  element  thickness.  This  result  was  not  used  because  the  thickness  of 
each  layer  was  assumed  to  be  significant  to  the  sensor  performance.  Decreasing  the 
thickness  would  have  altered  the  distance  between  the  various  parts  of  the  antenna, 
essentially  producing  a  different  antenna  element.  Also,  the  intent  was  to  validate 
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Figure  4.1  First  buckling  mode  shape  for  all  edges  simply  supported 


Figure  4.2  First  buckling  mode  shape  with  chordwise  edges  clamped  and  spanwise 
edges  simply  supported 

the  stiffness  of  the  current  model,  and  not  to  reproduce  the  exact  results  achieved 
by  Northrop  Grumman. 

4-2  Airfoil  Curvature 

The  radius  of  curvature  calculations  found  that  40%  of  the  bottom  surface 
of  the  wing  had  a  radius  of  curvature  greater  than  10  times  the  wavelength,  A. 
Approximately  71%  of  the  top  surface  of  the  wing  had  a  radius  of  curvature  that 
did  not  meet  the  length  criteria,  which  suggested  that  an  array  could  be  considered 
flat  only  if  the  local  radius  of  curvature  was  greater  than  10A.  A  series  of  plots  is 
shown  in  Figure  4.4.  The  top  plot  shows  the  shape  of  the  airfoil  as  a  reference,  the 
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Figure  4.3  Half  sine  wave  mode  shape 


middle  plot  is  the  absolute  value  of  the  radius  of  curvature  of  the  top  surface,  and  the 
bottom  plot  is  the  absolute  value  of  the  radius  of  curvature  of  the  bottom  surface. 
The  main  array  structure  was  assumed  to  fit  between  the  leading  and  trailing  edge 
spars.  Looking  at  the  plot  of  the  airfoil,  the  leading  edge  spar  was  at  an  x-coordinate 
of  11.96  inches,  or  12.2%  of  the  chord,  and  the  trailing  edge  spar  was  at  73  inches,  or 
74.5%  chord.  With  the  added  wing  deformation,  neither  the  top  nor  bottom  surface 
remains  planar  regardless  of  the  chordwise  curvature.  For  these  reasons,  the  flat 
array  approximation  was  not  considered  for  this  study. 

4-3  Pattern  Computation  with  NEC-Win  Plus+™ 

NEC-Win  Plus+™  was  also  used  to  obtain  radiation  patterns  for  the  confor¬ 
mal  array  of  dipoles  on  the  undeformed  wing  and  in  the  case  where  a  gust  load  was 
applied  to  the  wing.  NEC-Win  Plus-1- ™  will  only  plot  the  gain  pattern  of  the  array 
in  dB,  but  the  basic  shape  of  the  radiated  pattern  should  be  similar  whether  the 
plot  is  gain  or  simply  the  field  magnitude  plotted  in  W  The  radiation  pattern  of  a 
10x40  element  array  located  close  to  the  fuselage  was  calculated  for  the  undeformed 
wing  and  the  gust  load  case.  The  results  from  the  basic  array  analysis  are  shown  in 
Figures  4.5,  4.6,  4.7,  and  4.8. 
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Radius  (in)  Radius  (in)  Y  Coordinate  (in) 


Figure  4.4  Airfoil  (top),  top  surface  radius  of  curvature  (middle),  and  bottom 
surface  radius  of  curvature  (bottom) 
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Figure  4.5  Magnitude  of  <j)  component  of  E-field  radiation  in  from  10x40  array 
on  undeformed  wing  with  end-fire  phasing:  elevation  variation  (0  = 
180°) 


Figure  4.6  Magnitude  of  <j)  component  of  E-£eld  radiation  in  [^]  from  10x40  array 
on  undeformed  wing  with  end-fire  phasing:  azimuth  variation  (6  =  90°) 

These  plots  were  generated  in  MATLAB  showed  the  linear  magnitude  of  the  E-field 
for  variations  in  the  elevation  angle,  6 ,  and  azimuth  angle,  q i.  These  values  were 
compared  to  the  NEC  plots  of  power  gain  since  the  main  beam  pointing  angles 
should  have  been  the  same.  The  patterns  from  the  NEC-Win  Plns+™  analysis  are 
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Figure  4.7  Magnitude  of  0  component  of  E-field  radiation  in  from  10x40  array 
for  gust  load  case  with  end-fire  phasing:  elevation  variation  (0  =  180°) 


Figure  4.8  Magnitude  of  0  component  of  E-field  radiation  in  [h-]  from  10x40  array 
for  gust  load  case  with  end-fire  phasing:  azimuth  variation  (6  =  90°) 

shown  in  Figures  4.9,  4.10,  4.11,  and  4.12.  The  main  beam  magnitude  and  angular 
location  are  shown  in  the  ”  Pattern  Analysis”  box  in  each  figure. 

The  NEC-Win  Plus-1- ™  plots  for  elevation  variation  were  flipped  relative  to 
the  plots  generated  by  the  MATLAB  code.  The  main  beam  in  both  of  these  plots 
was  still  pointing  toward  +90°,  but  the  perspective  changed.  The  MATLAB  code 
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+-180 

Figure  4.9  NEC-Win  Plus+™  plot  of  gain  in  dB  for  the  undeformed  wing  case: 
elevation  variation  (0  =  180°) 


180 


Figure  4.10  NEC-Win  Plns+™  plot  of  gain  in  dB  for  the  undeformed  wing  case: 
azimuth  variation  ( 6  =  90°) 

viewpoint  is  essentially  the  root  of  the  wing  with  the  beam  shooting  off  to  the  left,  but 
the  NEC  viewpoint  is  from  the  direction  of  the  tip  of  the  wing.  Also,  the  direction  of 
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+-180 

Figure  4.11  NEC- Win  Plus-1- ™  plot  of  gain  in  units  of  dB  for  the  gust  load  de¬ 
formation  case:  elevation  variation  (0  =  180°) 


180 


Figure  4.12  NEC- Win  Plus-1- ™  plot  of  gain  in  units  of  dB  for  the  gust  load  de¬ 
formation  case:  azimuth  variation  (6  =  90°) 

180°  in  the  NEC  azimuth  plots  is  downward  instead  of  to  the  left  as  is  was  depicted 
in  the  MATLAB  plots. 
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Table  4.1  Comparison  of  NEC-Win  Plus+™  and  basic  array  theory  results  for 
angular  location  and  magnitude  of  E-field  maxima  for  undeformed  and 
gust  load  cases 


TECHNIQUE 

LOAD 

VARIATION 

MAX.  ANGLE 

MAG.  (dB) 

Array  Theory 

Undef  Wing 

e 

91.1° 

15.301 

NEC-Win  Plus+™ 

Undef  Wing 

9 

89.0° 

20.31 

Array  Theory 

Undef  Wing 

<t> 

179.91° 

15.30 

NEC-Win  Plus+riW 

Undef  Wing 

180.0° 

20.27 

Array  Theory 

Gust  Load 

e 

85.37° 

17.062 

NEC-Win  Plus+™ 

Gust  Load 

9 

92.0° 

22.26 

Array  Theory 

Gust  Load 

<t> 

181.05° 

15.29 

NEC-Win  Plus+™ 

Gust  Load 

</> 

181.0° 

22.21 

The  patterns  generated  by  NEC-Win  Plus-1- ™  were  more  refined  than  the 
plots  generated  using  basic  array  theory.  The  NEC  plots  show  deep  nulls  separating 
many  grating  lobes,  but  these  distinctions  are  lost  when  simple  held  summations 
are  used.  The  overall  patterns  show  some  similarities,  but  the  array  theory  patterns 
were  wider  than  the  NEC  patterns.  A  comparison  of  the  results  generated  by  the  two 
methods  are  listed  in  Table  4.3.  These  results  were  generated  for  a  10x40  element 
array  positioned  with  the  inboard,  chordwise  edge  at  the  root  of  the  front  wing. 
This  was  the  position  of  the  array  when  the  NEC  analysis  was  performed.  One 
run  in  NEC-Win  Plus+™  took  22.5  hours  for  a  10x40  element  array  on  a  1.7  GHz 
Windows  PC  with  256  MB  memory.  Also,  the  cases  that  were  run  validated  that  the 
MATLAB  code  was  a  decent  estimation  of  the  numerical  solution  for  the  undeformed 
wing,  but  it  was  less  reliable  for  predicting  the  main  beam  pointing  angle  as  the  wing 
deflection  increased. 

The  elevation  pointing  angle  calculated  by  NEC-Win  Plus+™  was  pitched 
up  2.1°  compared  to  the  array  theory  angle  for  the  undeformed  wing  case,  and  the 
magnitude  was  5  dB  higher.  In  the  azimuth  plane  the  angles  were  nearly  identical, 
but  the  magnitudes  differed  again  by  5  dB.  For  the  gust  load  case,  the  comparison  in 
elevation  angles  was  opposite  the  results  of  the  undeformed  wing.  The  array  theory 
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determined  that  the  main  beam  pitched  up  almost  6°  compared  to  the  undeformed 
wing,  but  NEC-Win  Plus+™  found  that  the  beam  was  actually  steered  downward 
3°  from  its  undeformed  wing  result.  The  azimuth  angles  were  essentially  the  same 
once  again  and  the  magnitude  differences  were  approximately  5  dB. 

4-4  Antenna  Results 

For  the  initial  antenna  analysis,  coordinates  for  the  deformed  array  extending 
40  elements  in  the  spanwise  direction  were  derived  from  a  2.5g  steady  maneuver 
load  case.  The  array  was  located  just  inboard  of  the  joint  section,  as  was  shown  in 
Figure  3.2.  This  wing  loading  produced  a  three-meter  vertical  deflection  at  the  tip 
of  the  wing,  but  only  a  one-meter  deflection  at  the  outboard  edge  of  the  antenna 
array.  When  the  load  case  was  changed  to  include  a  non-linear  deflection  for  a  gust 
load,  where  the  structure  was  sized  for  linear  analysis,  the  outboard  edge  of  the  array 
deflected  over  six  meters.  The  array  size  was  also  increased  to  60  elements  running 
in  the  spanwise  direction  to  see  if  a  larger  array  was  more  sensitive  to  the  same  load 
configurations.  The  plots  shown  in  this  section  are  the  linear  azimuth  and  elevation 
plots  of  the  E-field  horizontal  component  magnitude  for  both  of  these  load  cases  and 
for  the  undeformed  wing  case. 

Figures  4.13  and  4.14  are  the  undeformed  wing  array  radiation  patterns  with 
no  phasing  applied.  As  seen  in  Figure  4.13  for  the  undeformed  wing,  curvature  of 
the  antenna  as  it  conforms  to  the  airfoil  shape  in  the  chordwise  direction  causes  a 
shift  in  the  pointing  direction  of  the  main  beam  as  compared  to  a  planar  array.  If  the 
airfoil  was  flattened  (z-coordinates  set  to  zero),  the  theta  variation  pattern  looked 
like  Figure  4.15.  For  the  curved  array,  the  main  beam  shifts  clockwise  approximately 
seven  degrees. 

The  large  grating  lobe  that  appears  opposite  the  main  beam  in  these  figures  was 
expected  unless  phasing  was  applied  that  steered  the  main  beam  to  end-fire,  and  the 
interelement  spacing  was  considerably  less  than  half  a  wavelength.  It  was  discovered 
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Figure  4.13  Magnitude  of  (p  component  of  E-ficld  radiation  in  from  5x40  array 
on  undeformed  wing  with  zero  phasing:  elevation  variation  (0  =  180°) 


Figure  4.14  Magnitude  of  0  component  of  E-ficld  radiation  in  from  5x40  array 
on  undeformed  wing  with  zero  phasing:  azimuth  variation  (9  =  90°) 

that  simply  changing  the  phasing  to  an  end-fire  configuration  did  not  eliminate  the 
grating  lobe.  Therefore  the  number  of  dipoles  used  in  generating  these  plots  was 
increased  from  5  to  10  for  the  same  array  chord  length.  This  was  done  because  the 
spacing  for  five  dipoles  was  approximately  half  a  wavelength.  The  actual,  conformal 
elements  are  intended  to  have  only  one  main  beam  and  no  large  grating  lobes.  The 
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Figure  4.15  Magnitude  of  0  component  of  E-ficld  radiation  in  from  5x40  array 
on  flattened  airfoil  with  zero  phasing:  elevation  variation  (0  =  180°) 

CLAS  elements  accomplish  this  due  to  their  geometry,  but  very  thin  dipoles  at 
the  same  center  point  locations  do  not.  Increasing  the  number  of  dipoles  to  10 
corresponds  to  roughly  a  quarter  wavelength  spacing  for  a  frequency  of  430  MHz. 
Figures  4.16,  4.17,  and  4.18  show  the  undeformed  array  Great  Circle  pattern  cuts 
of  a  10x40  element  array  with  phasing  applied  to  steer  the  beam  to  an  end-fire 
configuration  and  a  smaller  interelement  spacing.  The  extremely  small  magnitude  of 
the  radiated  field  in  Figure  4.18  was  a  common  to  all  the  load  cases,  and  therefore 
only  one  plot  is  presented  for  this  pattern  cut,  which  pertains  to  an  elevation,  6 , 
variation  with  0  =  90°. 

The  effects  of  typical  wing  deformation  on  the  10x40  array  pointing  accuracy 
proved  to  be  significant  for  the  2.5g  steady  load  case,  and  even  more  so  for  the  gust 
load  case.  The  radiation  patterns  for  this  case  are  shown  in  Figures  4.19,  4.20,  4.21, 
and  4.22. 

While  the  differences  are  hard  to  see  on  these  plots,  the  angular  location  of  the 
maxima  changed  considerably.  The  undeformed,  2.5g  load,  and  gust  load  radia¬ 
tion  patterns  for  elevation,  6 ,  and  azimuth,  0,  variations  are  shown  in  Table  4.4, 
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Figure  4.16  Magnitude  of  <f>  component  of  E-field  radiation  in  [^]  from  10x40 
array  on  undeformed  wing  with  end-fire  phasing:  elevation  variation 

{<j>  =  180°) 


Figure  4.17  Magnitude  of  <fi  component  of  E-field  radiation  in  [^]  from  10x40 
array  on  undeformed  wing  with  end-fire  phasing:  azimuth  variation 
(i 9  =  90°) 

which  lists  the  angular  location  of  the  maxima,  linear  magnitudes  in  units  of  and 
magnitudes  in  dB. 

Considering  the  undeformed  wing  case  as  the  baseline,  the  maximum  value  of 
the  E-field  magnitude  for  the  elevation  plots  remained  virtually  the  same  when  the 
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Figure  4.18  Magnitude  of  0  component  of  E-field  radiation  in  [bj  from  10x40 
array  on  undeformed  wing  with  end-fire  phasing:  elevation  variation 
(0  =  90°) 


0 


Figure  4.19  Magnitude  of  0  component  of  E-field  radiation  in  [b]  from  10x40 
array  for  2.5g  load  case  with  end-fire  phasing:  elevation  variation  (0  = 
180°) 

2.5g  load  was  applied,  but  it  decreased  by  0.031  b  for  the  gust  load  case.  In  the 
azimuth  plots,  the  maximum  value  increased  by  0.002  b  for  the  2.5g  load,  but  the 
magnitude  for  the  gust  load  decreased  by  0.159  b.  The  angular  location  of  the  E- 
ficld  elevation  maximum  magnitude  decreased  for  both  cases  of  applied  loads.  The 
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Figure  4.20  Magnitude  of  4>  component  of  E-field  radiation  in  from  10x40 
array  for  2.5g  load  case  with  end-fire  phasing:  azimuth  variation  (6  = 
90°) 


Figure  4.21  Magnitude  of  <fi  component  of  E-field  radiation  in  [^]  from  10x40 
array  for  gust  load  case  with  end-fire  phasing:  elevation  variation 
(0  =  180°) 

gust  load  case  was  the  most  severe  of  the  two  loads,  and  caused  the  main  beam  to 
shift  up  9.17°.  In  the  azimuth  plots,  the  2.5g  load  caused  no  shift  in  the  pointing 
angle,  but  the  gust  load  caused  the  main  beam  to  shift  to  the  left  1.72°. 
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Figure  4.22  Magnitude  of  <f>  component  of  E-field  radiation  in  from  10x40 
array  for  gust  load  case  with  end-fire  phasing:  azimuth  variation  (6  = 
90°) 


Table  4.2  Angular  location  and  magnitude  of  front  wing  10x40  element  array  E- 
field  maxima  for  undeformed,  2.5g  load,  and  gust  load  cases 


LOAD  CASE 

VARIATION 

MAX.  ANGLE 

mag.  m 

V  777,  / 

MAG.  (dB) 

Undeformed  Wing 

e 

92.25° 

33.890 

15.301 

2.5g  Maneuver  Load 

9 

88.55° 

33.888 

15.30 

Gust  Load 

9 

81.80° 

33.857 

15.297 

Undeformed  Wing 

(t> 

179.91° 

33.880 

15.30 

2.5g  Maneuver  Load 

</> 

179.91° 

33.882 

15.30 

Gust  Load 

0 

181.63° 

32.723 

15.279 

The  array  was  lengthened  to  60  elements  in  the  spanwise  direction  to  see  if 
a  longer  array  would  have  more  noticeable  degradation  in  pointing  accuracy  due 
to  wing  deformation.  The  radiation  pattern  cuts  from  the  10x60  element  array  on 
an  undeformed  wing  are  shown  in  Figures  4.23  and  4.24.  The  results  are  shown  in 
Table  4.4. 

Again  considering  the  undeformed  wing  case  as  the  baseline,  the  maximum 
value  of  the  E-field  magnitude  for  the  elevation  plots  decreased  by  only  0.006  ^ 
when  the  2.5g  load  was  applied,  but  it  decreased  by  0.085  —  for  the  gust  load  case. 
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Figure  4.23  Magnitude  of  <f>  component  of  E-field  radiation  in  [^]  from  10x60 
array  on  undeformed  wing  with  end-fire  phasing:  elevation  variation 

{<j>  =  180°) 


Figure  4.24  Magnitude  of  <fi  component  of  E-field  radiation  in  [^]  from  10x60 
array  on  undeformed  wing  with  end-fire  phasing:  azimuth  variation 
(i 9  =  90°) 

In  the  azimuth  plots,  the  maximum  value  decreased  by  0.025  ^  for  the  2.5g  load, 
and  the  magnitude  for  the  gust  load  decreased  by  0.286  The  angular  location 
of  the  E-field  elevation  maximum  magnitude  once  again  decreased  for  both  cases  of 
applied  loads.  The  gust  load  case  caused  the  main  beam  to  shift  up  8.02°.  For  the 
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Table  4.3  Angular  location  and  magnitude  of  front  wing  10x60  element  array  E- 
field  maxima  for  undeformed,  2.5g  load,  and  gust  load  cases 


LOAD  CASE 

VARIATION 

MAX.  ANGLE 

mag.  m 

\  m  / 

MAG.  (dB) 

Undeformed  Wing 

e 

90.53° 

50.839 

17.062 

2.5g  Maneuver  Load 

e 

87.40° 

50.833 

17.061 

Gust  Load 

9 

81.23° 

50.754 

17.055 

Undeformed  Wing 

<t> 

179.91° 

50.837 

17.062 

2.5g  Maneuver  Load 

<t> 

180.48° 

50.812 

17.060 

Gust  Load 

182.20° 

50.551 

17.037 

azimuth  plots,  the  2.5g  load  caused  a  0.57°  shift  in  the  pointing  angle,  and  the  gust 
load  caused  the  main  beam  to  shift  to  the  left  2.29°. 

A  10x60  element  array  was  placed  on  the  aft  wing.  The  radiation  pattern  cuts 
from  the  array  on  an  undeformed  wing  are  shown  in  Figures  4.25  and  4.26.  The 
patterns  due  to  the  2.5g  maneuver  load  are  shown  in  Figures  4.27  and  4.28,  and  for 
the  gust  load  the  linear  E-field  plots  are  shown  in  Figures  4.29  and  4.30  The  results 
are  shown  in  Table  4.4. 


0 


Figure  4.25  Magnitude  of  (f>  component  of  E-held  radiation  in  from  10x60  ar¬ 
ray  on  undeformed  aft  wing  with  end-fire  phasing:  elevation  variation 
(0  =  180°) 


Deviation  from  the  undeformed  wing  case  showed  that  the  maximum  value  of 
the  E-held  magnitude  for  the  elevation  plots  decreased  by  only  0.005  ^  when  the 
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Figure  4.26  Magnitude  of  <f>  component  of  E-field  radiation  in  from  10x60 
array  on  undeformed  aft  wing  with  end-fire  phasing:  azimuth  variation 
(< 9  =  90°) 


Figure  4.27  Magnitude  of  <fi  component  of  E-field  radiation  in  [^]  from  10x60 
array  on  aft  wing  for  2.5g  load  case  with  end-fire  phasing:  elevation 
variation  (0  =  180°) 

2.5g  load  was  applied,  and  it  decreased  by  0.009  ^  for  the  gust  load  case.  So  there 
was  essentially  no  change  in  magnitude.  In  the  azimuth  plots,  the  maximum  value 
decreased  by  0.035  —  for  the  2.5g  load,  and  the  magnitude  for  the  gust  load  decreased 
by  0.146  — .  The  angular  location  of  the  E-held  elevation  maximum  magnitude 
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Figure  4.28  Magnitude  of  0  component  of  E-field  radiation  in  from  10x60 
array  on  aft  wing  for  2.5g  load  case  with  end-fire  phasing:  azimuth 
variation  ( 6  =  90°) 


Figure  4.29  Magnitude  of  0  component  of  E-field  radiation  in  [^]  from  10x60 
array  on  aft  wing  for  gust  load  case  with  end-fire  phasing:  elevation 
variation  (0  =  180°) 

decreased  for  both  cases  of  applied  loads.  The  gust  load  case  caused  the  main  beam 
to  shift  up  9.14°.  For  the  azimuth  plots,  the  pointing  angle  remained  virtually 
unchanged. 
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Figure  4.30  Magnitude  of  0  component  of  E-field  radiation  in  from  10x60 
array  on  aft  wing  for  gust  load  case  with  end-fire  phasing:  azimuth 
variation  ( 9  =  90°) 

Table  4.4  Angular  location  and  magnitude  of  aft  wing  10x60  element  array  E-field 


maxima  for  undeformed,  2.5g  load,  and  gust  load  cases 


LOAD  CASE 

VARIATION 

MAX.  ANGLE 

mag.  m 

V  m.  / 

MAG.  (dB) 

Undeformed  Wing 

e 

271.0° 

50.851 

17.063 

2.5g  Maneuver  Load 

9 

273.87° 

50.846 

17.063 

Gust  Load 

9 

261.84° 

50.842 

17.061 

Undeformed  Wing 

0 

O 

O 

o 

50.847 

17.063 

2.5g  Maneuver  Load 

0 

O 

o 

o 

50.812 

17.060 

Gust  Load 

0 

359.82° 

50.701 

17.050 

Increasing  the  size  of  the  array  caused  the  magnitude  of  the  electromagnetic 
field  to  be  greater,  but  the  magnitude  was  slightly  more  sensitive  to  deformation 
than  it  was  for  the  shorter  array.  Comparing  the  pointing  accuracy  of  the  two  arrays 
showed  that  the  larger  array  had  a  slight  advantage  in  the  elevation  variation,  but 
its  main  beam  experienced  a  larger  shift  in  the  azimuth  plane.  Neither  array  was 
able  to  maintain  the  desired  end-fire  pointing  configuration  of  9  =  90°  and  0  =  180°. 
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4-5  Aileron  Effectiveness 

The  control  surface  was  found  to  be  only  5%  effective  in  generating  lift  to  trim 
the  aircraft  in  pitch.  This  confirmed  that  an  inboard  control  surface  or  aft  wing  twist 
ought  to  be  used  to  trim  this  aircraft  for  level  amneuvers.  The  roll  effectiveness  of 
the  aileron  was  also  evaluated  for  various  Mach  numbers.  At  the  intended  loiter 
Mach  number  of  0.6  (V  =  204.17  — ),  the  non-dimensional  roll  rate  was  found  to  be 
slightly  negative  (p  =  — 5.01a;10”4),  which  meant  that  aileron  reversal  had  occurred. 
The  dynamic  pressure  was  varied  slightly  to  bracket  the  reversal  speed.  At  a  velocity 
of  203.64  —  the  non-dimensional  roll  rate  was  positive  (p  =  3.31a;10~4),  which  meant 
that  the  aileron  was  rolling  the  aircraft  in  the  intended  direction. 
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V.  Conclusion 


This  study  shows  that  the  radiation  patterns  of  an  arrays  that  conform  to  the  surface 
of  the  front  and  aft  wing  sections  of  the  joined  wing  undergo  significant  distortion  dne 
to  typical  wing  deflections.  For  a  highly  deformed  wing  dne  to  a  gust  load,  simple 
phasing  of  the  array  was  not  adequate  to  keep  the  beam  pointed  in  the  desired 
direction.  Even  for  a  steady  2.5g  maneuver  load,  the  change  in  the  main  beam 
elevation  angle  was  substantial.  Active  control  by  sensing  wing  deformations  may 
help  compensate  for  some  of  the  error  in  beam  steering,  but  structural  modifications 
may  also  be  necessary  for  the  configuration  described  in  this  effort. 

Comparison  of  the  patterns  with  a  Method  of  Moments  solution  showed  a  6.6° 
degree  difference  in  the  main  beam  elevation  angle  for  a  gust  load,  but  only  a  2° 
difference  for  the  undeformed  wing  case.  The  azimuth  angles  for  all  of  the  load 
cases  were  essentially  the  same.  The  MATLAB  code  can  be  easily  incorporated  into 
an  aircraft  design  environment  where  a  structural  engineer  can  use  this  information 
to  design  aircraft  around  a  large  sensor  suite  such  as  the  one  presented  in  this 
effort.  The  NEC-Win  Plus+™  software  would  be  impossible  to  incorporate  into  the 
AVT1E  environment  because  each  source  has  to  be  placed  and  configured  manually. 
Therefore,  the  method  presented  provides  a  quick  and  reasonable  estimation  of  the 
conformal  array  radiation  pattern  of  an  array  placed  anywhere  on  the  undeformed 
wing,  but  it  is  less  reliable  as  the  wing  deflection  increases. 

A  control  surface  on  the  wing  tip  section  of  this  joined  wing  was  not  effective 
for  rolling  the  aircraft,  because  roll  reversal  occurs  just  before  Mach  0.6  at  sea  level 
conditions.  The  aileron  was  also  not  effective  for  pitching  the  aircraft.  Therefore, 
aft  wing  twist  was  justified  to  trim  the  aircraft  for  pitch. 
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5.1  Future  Research 


There  are  several  areas  of  this  research  that  were  either  not  attempted  due  to 
time  constraints  or  need  further  development.  These  topics  are  listed  below: 

1.  Increase  the  fidelity  of  the  current  finite  element  composite  model  for  a  more 
exact  design. 

2.  Look  at  a  vertical  dipole  or  some  other  element  (square-loop,  patch,  etc.). 

3.  Model  the  actual  antenna  elements  or  incorporate  contractor  data  into  the 
methodology  presented  here. 

4.  Place  sensor  arrays  on  the  bottom  surfaces  of  both  the  front  and  aft  wing 
sections  inboard  of  the  joint,  and  couple  top  and  bottom  arrays  together 

5.  Look  at  a  smaller,  high  frequency  antenna  out  on  the  tip  section  of  the  wing. 

6.  Optimize  the  various  configurations  of  the  antenna  arrays  and  control  surfaces. 
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Appendix  A.  LRN-1015  Data  Points 


UPPER 

LOWER 

X 

Z 

X 

Z 

0 

0 

0 

0 

0.001621 

0.01707 

0.001621 

-0.001976 

0.006475 

0.026197 

0.006475 

-0.00493 

0.014529 

0.036312 

0.014529 

-0.007513 

0.025732 

0.04699 

0.025732 

-0.01001 

0.04001 

0.057142 

0.04001 

-0.012525 

0.057272 

0.066677 

0.057272 

-0.014983 

0.077405 

0.075746 

0.077405 

-0.017226 

0.100279 

0.084202 

0.100279 

-0.019318 

0.125745 

0.092004 

0.125745 

-0.021214 

0.153638 

0.099088 

0.153638 

-0.022877 

0.183777 

0.105433 

0.183777 

-0.024304 

0.215968 

0.110974 

0.215968 

-0.025466 

0.25 

0.115663 

0.25 

-0.026357 

0.285654 

0.119458 

0.285654 

-0.026977 

0.322698 

0.122299 

0.322698 

-0.027302 

0.360891 

0.124114 

0.360891 

-0.02733 

0.399987 

0.124809 

0.399987 

-0.027065 

0.439732 

0.124264 

0.439732 

-0.026505 

0.479867 

0.122351 

0.479867 

-0.025652 

continued  on  next  page  . . . 
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0.520133 

0.118919 

0.520133 

-0.024494 

0.560268 

0.113813 

0.560268 

-0.023028 

0.600013 

0.106947 

0.600013 

-0.021274 

0.639109 

0.098407 

0.639109 

-0.019239 

0.677302 

0.088455 

0.677302 

-0.016865 

0.714346 

0.077492 

0.714346 

-0.014081 

0.75 

0.06623 

0.75 

-0.010938 

0.784032 

0.055549 

0.784032 

-0.007663 

0.816223 

0.046102 

0.816223 

-0.004646 

0.846362 

0.037889 

0.846362 

-0.00213 

0.874255 

0.030724 

0.874255 

-0.000215 

0.899721 

0.024801 

0.899721 

0.001069 

0.922595 

0.020019 

0.922595 

0.001761 

0.942728 

0.015794 

0.942728 

0.001957 

0.95999 

0.01182 

0.95999 

0.001792 

0.974268 

0.008107 

0.974268 

0.001378 

0.985471 

0.004834 

0.985471 

0.000884 

0.993525 

0.002244 

0.993525 

0.000429 

0.998379 

0.000577 

0.998379 

0.000113 

1 

0 

1 

0 
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Appendix  B.  THE  FOURIER  TRANSFORM 

This  information  was  adapted  from  Dr.  Andrew  Terzuoli’s  Antennas  I  class  notes. 


Take  two  arbitrary  1-D  spaces  W  and  S,  where  W  and  S  are  real  variables  such  that 
their  product  (WS)  is  unitless.  These  spaces  are  called  inverse  spaces,  with  complex 
functions  f(W),  F_(S)  that  can  be  related  by  the  Fourier  Transform  given  by: 


f(W) 

F(S) 


F(S )  ej2n(WS)dS 
f_(W)e-j2n{ws)dW 


Now  apply  to  radiation  and  scattering  by  a  current,  where 


1r  (9,  1r  (r,  u)  ePLW+VWd \T 


\?\ 


and 


AR(rs,u )  =  /f(w)I(r,w)Jfl(M,w). 


In  these  equations,  7 
of  the  radiation,  and 


a  +  j(3,  where  a  is  the  attenuation  and  (3  is  the  propagation 


T  = 


e~T 
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An  engineering  approximation  is  made  here  that  ignores  the  transverse  attenuation 
in  the  integral.  Therefore,  in  the  integral  there  is  only  transverse  resonance  and 
7  =  jfi.  Then 

Ir'l 

where  j/3(u))rl  ■  r(6,  R)  is  a  real  number. 
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The  source  is  now  restricted  to  rectangular  variables,  thus  let 

rR  —  x  x  +  y  y  +  z  z 

and  the  dot  product,  r'R  ■  r,  becomes 

r'n  ■  r  =  x'(x  ■  r)  +  y\y  ■  r)  +  z\z  ■  f) 

where  the  dot  products  of  this  equation  are  the  direction  cosines: 

(x  ■  r )  =  cos  9X 
(y  ■  r)  =  cos  9y 
(z  ■  r)  =  cos  9 z. 

If  now  transformed  to  spherical  coordinates,  the  direction  cosines  become 

cos  9 x  =  sin  9  cos  0 
cos  9y  =  sin  9  sin  cj) 
cos  9  z  =  cos  9. 

Then  defining 

(3X  =  (3  cos  9X 
(3y  =  f3  COS  9y 
(3Z  =  (3  cos  9Z 
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/ 3  times  the  vector  dot  product  becomes 


PifR  ■  f)  =  x'(3x  +  y'py  +  z'pz  =  rR  ■  0,  u). 


This  real,  unitless  dot  product  is  substituted  back  into  the  integral  to  give 


Jr  {0,4>,u) 


which  is  a  three-dimensional  Fourier  Transform. 
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